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1.0 SUMMARY 


During recent years there has been a growing interest in new air- 
craft propulsion systems such as powered-lift devices and high- 
speed propellers (propfans) for both commercial and military 
applications. Studies of the systems have indicated that they 
will generate high noise levels, particularly at low frequencies, 
in the fuselage interiors. At the same time there has been an 
increasing demand on the National Aeronautics and Space Adminis- 
tration to provide technical support to the general aviation 
industry, and one problem identified in this field is the reduc- 
tion of cabin noise levels. In an attempt to solve some of the 
complicated problems involved in the control of airplane interior 
noise, NASA Langley Research Center embarked on an extensive 
program to identify the important parameters associated with 
acoustic transmission (airborne and structureborne ) into airplane 
interiors and to determine appropriate noise control methods 
which are acceptable with regards to weight and space. As a part 
of this program, an aircraft interior noise prediction model is 
being developed by Bolt Beranek and Newman Inc. (BBN), the 
purpose being to undertake the required analyses and integrate 
the technologies and information needed to understand sound 
transmission through a fuselage wall into an aircraft cabin, both 
for future development studies and for the immediate need of 
accurate prediction of interior levels (more properly stated, 
for prediction of the noise isolation characteristics of the 
fuselage structures). The effort has been divided into three 
phases. These are: 



Phase I. 


Analytical modeling of the transmission problem 
(interaction of the structures with the exterior 
and interior acoustic fields) and preliminary 
verification studies using, as a test article, an 
unpressurized unstiffened cylinder. 

Phase II. 

Development of the general aircraft interior noise 
model and basic master computer program with appli- 
cation to advanced test articles (for instance, a 
stiffened cylinder with interior partitions and 
simulated trim) . 

Phase III. 

Completion of the analytical models and software 
with application to an actual aircraft fuselage, 
verification tests, comparisons, refinements, and 
documentation of the finalized model and software. 

This report presents the results of the Phase I studies. The 
basic theoretical work required to understand sound transmission 
into an enclosed space (that is, one closed by the transmitting 
structure) is developed for random pressure fields and for har- 
monic (tonal) excitation. The analysis is used to predict the 
noise reduction of the test article (i.e., the unpressurized 
unstiffened cylinder) and, also, the interior response of the 
cylinder given a tonal (plane wave) excitation. Comparisons of 
predictions and measurements are made and a detailed analysis of 
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the transmission is presented. This latter study is perhaps the 
most informative aspect of the exercise undertaken in Phase I. 

In addition to the above, results for tonal (harmonic) mechanical 
excitation are considered. 

The work was accomplished as a joint effort by BBN/Los Angeles 
and NASA Langley Research Center. L. D. Pope served as BBN 
program manager and William H. Mayes was LaRC technical monitor. 
The authors wish to thank Mr. Mayes and Conrad Willis of NASA 
LaRC for their contributions to the experimental program; also 
Allan G. Piersol of BBN who performed the statistical analysis 
of the noise reduction data, and John F. Wilby of BBN for general 
assistance . 
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2.0 INTRODUCTION 


The present study has the specific goal of developing an analy- 
tical model which can be used to predict the sound levels in a 
cylinder for several excitation inputs and interior configura- 
tions. However, this objective has to be viewed in terms of the 
overall objective of the program of which this present study 
forms only the first phase. The complete program has the objec- 
tive of developing an analytical model to describe the sound 
field inside an airplane cabin. Consequently, the analytical 
model resulting from the present study must not only predict the 
sound levels in the unstiffened and unpressurized test cylinder 
(which is an idealized fuselage structure), but it must be 
adaptable to the more complicated situation present in a realis- 
tic airplane fuselage. Only then can this study be viewed as a 
first phase in the larger program. 

The two main efforts of Phase I are the construction of the 
analytical models and the experimental validation of those 
models on an unstiffened cylinder. The work effort under 
Phase I and envisaged for Phases II and III of the present 
program can be considered in terms of the interior noise re- 
search and development program underway at Langley Research 
Center. This overall program is concentrating on two areas: 

(1) the accumulation of experimental data which can be incor- 
porated into a prediction methodology, and (2) development of 
the prediction methodology itself. The work has a bias towards 
general aviation aircraft. 

Interior noise levels have been measured by Catherines and 
Mayes [1] in single- and twin-engined propeller-driven light 
aircraft. These studies define to a large extent the nature of 
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the acoustic environment in light aircraft and establish the need 
for analytical prediction methods which will eventually be use- 
ful in the design of quiet interiors. Sound levels in the cabins 
of the aircraft were found to be dominated by propeller and 
engine exhaust noise, the spectra being essentially tonal in 
nature. In [2] the measured exterior spatial distributions of 
the acoustic pressures showed the dominance of propeller noise 
over the forward part of the fuselage, and exhaust noise over the 
aft part. Reduction of interior noise is clearly desirable since 
levels in light aircraft range from 85 to 105 dBA. Barton [ 3 ] 
developed a method for predicting interior noise in STOL aircraft 
and, utilizing a sidewall noise reduction model which likely is 
conservative, predicted levels of 90 - 105 dBA and 88 - 92 dBA 
in aircraft with USBF and EBF systems respectively. 

Howlett et al . , [4] measured the noise reduction of a light 
aircraft in a reverberant room. Using the Cockburn and Jolly 
model [51, a comparison of predicted and measured noise reduc- 
tion was made. Agreement was not very satisfactory. 

In addition to the above studies, a series of progressively more 
complicated analysis and data reduction investigations has been 
completed under LaRC sponsorship. These include studies by 
Koval [6, 7, 8], Vaicaitis [ 9 ] and Rennison and Wilby [10,11] 
among others. 

2 . 1 Candidate Approaches for Analytical Model 

Before developing an analytical model for the present study, it 
is useful to review the technology available and to consider 
its applicability to the Phase I studies, and to subsequent 
Phase II and Phase III work. In general, it can be stated that 
much of the technological base required for the prediction of 
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the test cylinder response and acoustic radiation is available 
already. Potential analytical techniques which could be used 
for the present predictions may be generally grouped into two 
categories, each of which may be broken down further, basically 
in terms of computational schemes. When making predictions, 
such as required here, with potentially high acoustic and struc- 
ture modal densities, it is important to keep in mind that the 
choice of computational scheme is significant in the overall 
problem solving procedure. The two types of analyses are: 

1. Basic Normal Mode Approach 

(Pure tone and noise formulations; simple structures and 
geometrically simple cavities are specified). 

(a) Computations performed by straightforward summation 
over acoustic and structure modes. 

(b) Calculations performed by selecting the most likely 
responsive modes and by applying a convergence 
criterion to limit the number of summations. 

2. Power Balance Approach 

(a) Statistical energy analysis (high frequency band average 

predictions for noise). 

(b) Generalized power flow formulation (harmonics or noise, 
complex structures and geometrically complex cavities). 

Each of the above techniques has been used at one time or another 
in studies similar to that of the present Phase I. 

The Type 1(a) analysis and computation scheme might best be 
described as the brute force method. In principle, the exact 
dynamics of the structure have been developed by Pope [12], 
Cockburn and Jolly [ 5 ] and Pope [12] have solved the finite closed 
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cylinder problem for the case of a reverberant field (and in 
general for any homogeneous stationary random pressure field). 

The computer software procedures used in those studies are 
limited because the computation time requirements can be large, 
since summations over large numbers of acoustic and structure 
modes can be required. 

The analytical method of Type 1(b) is identical to Type 1(a) and 
might be typified by work such as that of Plotkin et al., [13]- 
The significant difference in this approach and 1(a) is the 
calculation procedure, which reduces computational times. The 
work in [ 13 ] represents an advance in the computational scheme 
for high modal density systems. Nevertheless, in general, the 
normal mode approach alone leaves much to be desired because it 
does not create a substantial understanding of the transmission 
phenomenon. That is, it gives answers but does not provide much 
auxiliary physical information to discern why certain results 
are computed. 

Type 2 analyses are significantly different from Type 1 in that 
net inflowing acoustic power is calculated and equated to net 
internal dissipation to obtain the interior levels. Early 
applications of this approach utilized statistical energy analy- 
sis, the Type 2(a) procedure. Analyses by Lyon [14], White [15], 
Scharton [16], and Eichler [17] have great utility but they all 
have the severe handicap that they are not applicable to low 
frequencies . 

Type 2(b) is a technical approach developed recently at BBN by 
-Pope and Wilby [18, 19, 20], for the acoustic environment in the 
payload bay of the space shuttle orbiter vehicle, that enables 
the use of power flow concepts at both low and high frequencies, 
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thereby overcoming a major disadvantage of the statistical energy 
analysis method. Essentially, beginning with modal analysis, a 
system of power flow and dissipation expressions has been 
developed which encompasses the normal mode approach at low 
frequencies and a modified form of statistical energy analysis 
at high frequencies. The technique considers power transmission 
by resonant and nonresonant structure modes into resonant and 
nonresonant acoustic modes. It can handle random and pure tone 
excitations at low or high frequencies, provide band-limited or 
harmonic predictions, and be converted into computer software 
which performs calculations over the entire frequency range of 
concern faster than can the method of [ 5 ]. It also provides 
auxiliary information, such as radiation efficiencies of struc- 
tures, band-averaged joint acceptances, etc., and can utilize 
modal data generated by different methods such as closed form 
equations or finite-element analyses. The method is not limited 
to simple structural models and it can handle complex receiving 
spaces. Moreover, its usefulness and accuracy have been demon- 
strated already in application to a structure as complex as the 
mid-fuselage of the space shuttle orbiter vehicle [21]. It is this 
latter technique which has been selected for the present study. 

2 . 2 Report Organization 

With these preliminaries behind, a brief overview of the organi- 
zation of this report may be useful. Basically this report 
presents results of (1) analytical derivations, (2) experimental 
tests, and (3) verification studies (comparisons). 

In the analytical development, there are always two fundamental 
requirements. The models have to have the capacity for applica- 
tion to a realistic fuselage, and, for the immediate needs of the 
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Phase I study, they must be specialized for application to the 
test article. Because of the first requirement, there is, of 
necessity, a certain abstractness in the fundamental developments. 
Every attempt has been made to allow the reader to follow the 
mathematics without excessive difficulty. Application of the 
general results to the case of the test cylinder will, hopefully, 
then reveal considerable insight into the meaning and usage of 
the models. It is important to understand this application 
because there will be a similarity involved when the models are 
used in reference to an actual fuselage. 

The parts of the report concerned with the presentation of the 
experimental work describe the purpose of the tests, data acqui- 
sition and reduction techniques, and give results usually as raw 
data, and in format as required for the comparison studies. 

Comparisons are of different types ranging from the "weak" quali- 
tative type such as visual comparisons of graphs of data versus 
predictions, to the "strong" statistical type of comparison deemed 
useful and necessary for verification purposes. 
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3.0 ANALYTICAL MODELS 


3 . 1 Generalized Power Flow Approach 

In the preceding section, a broad review of the general analysis 
techniques and calculation methods which are available to under- 
take Phase I, II, and III studies was given. In the review, the 
generalized power flow approach (Type 2(b)) was identified as the 
method most suitable to all three phases of the overall program. 
Thus it is appropriate to consider this approach in detail so 
that the component parameters can be identified. For purposes 
of discussion, it is assumed that the excitation is broadband 
acoustic noise, but the approach is not restricted to this 
excitation, as will be shown. 

The power balance method of analysis can be basically described 
with the equation 


W. = W, . 
in diss 


( 1 ) 


This equation simply states that the net time-averaged power, 
W^ n , flowing into the interior of the cylinder must equal the 
net time-averaged power, W diss , dissipated on the interior walls. 
By expressing the inflowing power in terms of the exterior ex- 
citing pressure and the dissipated power in terms of the inter- 
ior response pressure and equating the expressions according to 
Eq. (1), one can solve for the interior pressure. Thus, if 
these expressions are available, one can write down the solution 
to most sound transmission problems immediately. 


A primary feature of the generalized power flow approach adopted 
here is that integration of spectral components has been performed 
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analytically to achieve band-limited levels, i.e.. 


W 


in 


= / Re [ v O u) ] 

Am 


dm 


(2a) 


W 


diss 


/ Re [ W abs ( “ ) ] 
Am 


dm 


(2b) 


In the above expressions W ^^( w ) is the complex spectral density 
of the power radiated by the structure into the interior acoustic 
space and ^^(m) is the spectral density of power absorbed on 
the inner wall of the space from the interior acoustic field. 

Note here the similarity to statistical energy analysis (SEA has 
a band-limited power basis). However, Eqs. (2) have applicability 
at all frequencies and can be applied even when there are no 
acoustic or structure modes in the band Am. 


Power flow into each individual acoustic mode is computed, 

whether the mode is resonant in the band or not. If the mode is 

resonant in the band, one can speak of the "resonant response" of 

the mode; if the mode is resonant outside the band, the 

"nonresonant response" of the mode -in the band Am. Similarly 

power flowing out of the acoustic field to the inside wall is 

computed mode-by-mode. By equating the power into a mode to the 

power out, the space average mean square modal pressure for 

mode n, limited to band Am, is obtained, i.e., <p 2 > , . The band- 

n s j u 

limited, space-average mean-square pressure in the interior is the 
sum 



> 


s , t 


y <p 2 > 

Z, n s,t 


( 3 ) 
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Here n is the complete set of acoustic modes, i.e., those reso- 
nant inside Aw (denoted by the symbolism, neAw), as well as those 
resonant above and below Aco (n^Aw). At this point it should be 
noted that the letter "n" always denotes acoustic modes and the 
letter "r" structure modes, usually being used as subscripts, 
superscripts, or indexes. 

3.2 Solution of a General Sound Transmission Problem 


As previously stated, much of the technological base for the 
predictions of interest here has been developed [19,20]. Rather 
than repeat the derivations given in these two papers, the 
purpose of this section is to summarize and use the results. It 
is not absolutely necessary (but it will be helpful) for the 
reader to go through the papers before continuing with this 
section. The highlights of the developments are discussed 
below. 

The concern is with a closed volume V, into which sound energy 
is propagating. The excitation of the space occurs at its 
boundary with the enclosing structure which vibrates due to an 
exterior source. For the present, the excitation can be con- 
sidered to be a broadband random acoustic field. For the speci- 
fic case where the noise reduction is desired, the exterior field 
is taken to be diffuse (reverberant). 

Consider first the spectral density of the power that flows into 
the cavity. This quantity is obtained by integrating, over the 
transmitting area, the cross-power spectral density of interior 
pressure on the wall with the wall velocity at that position. 

wj^(u) = f Sp 1 v(x,ai)dx (4) 

A 
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The development of this expression into the desired form is 
presented in Sections II and III of Reference [19], resulting in 
Eq. (13) of Section III of [19], which is Eq . (5) below. 


r -t 2pto 3 A 4 S .-.(w) e A A 

Re[w ln j(io)l = Y n n n 

L rad J v 4 - I X 2 -k 2 1 2 


n 


n 


* 2 


f ' 2 ( n , r ) .j 2 ( a) ) 


:M r Y r -pa) 2 (J rr + I rr )| 2 


Nomenclature here is consistent with that in References 
This result is substituted into Eq . (2a), to obtain 


W. 

in 


pc 2 A 4 S , , ( to ) — » j 2 ( w ) 

o pbl Y _ r 


V 


„ M z u> 4 
r r r 


e n f ' 2 (n,r) 
V n n 

x z — ^ — 

^ n 


/ 

A u 


a) 3 da) 


{ [1- (co 2 /oj 2 ) ] 2 + n 2 H[l-(a) 2 A3 2 )] 


n 


rr 


Eq. (6) is the fundamental result for the input power, 
remains to obtain the desired left-hand side of Eq. (1) 
integrate the expression (Eq. (6) is the same as Eq . (3 
Section I of [20] ) . 


(5) 

[ 19 , 20 ]. 

(6) 

2 + njf) 

All that 
is to 
> of 
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Different forms for the inflowing power, , are obtained in 
References [19] and [20]; e.g., compare Eq . (19) of Section III 
of [19] to Eq . (5) of Section II of [20]. This difference has to 

do with how the integration is performed, i.e., with or without 
restrictions being placed upon the locations of the resonance 
frequencies of the modes (either inside or outside the band Aw). 
In Reference [20], the integration has been performed without 
restriction on mode locations. Thus Eq. (5) of Section II of 
[20] is the general result. The restriction njefAw (i.e., that w n 
lie outside Aw) can be lifted since it is artificial. Therefore, 
the required form of the left-hand side of Eq. (1) is 


W. = 
in 


< Pb3 > 

c w 
w 


2ttA' 

~w 


( mA V V 

XT/ t ml 2 


n 


e n 
n n 

r 


j 2 (u>)f ' 2 (n,r) i2 w 2 w 2 


M 2 

r 


D 


nr 



'2c (b -b )-b (c -c ) 

m TO Vv ' vn V vi ' 


♦I n ' r — ~ n " r — ) arctan ln 

Hn / n V H ) r 


( 


2c (b -b )-b (c -c ) 

VJ -VO ' v> V TO m ' 


j., rnr r ' n r ' , 

+ 1 1 arctan 


Hr] ■■ 2 


w 

r r 


) arctan l 

r J) 


(7) 


where 


D nr = (c r" c n )2 + (b n“ b r } (b n C r"Vn ) 


with 


b = - 
n 


2w , b 
n r 


-2w 2 , 

r * 


“ + hM , c = w“* (1 + t)D 


n n 


n 
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and for n or r = j , 


/ | ( 1+c /2 ) 4 co 4 +b . ( 1 + c /2) 2 co 2 +c.| \ 
l n _ 1n a .1 QJ J 1 I 

j \ (1-c /2 ) 4 co 4 +b^ ( 1-c /2) 2 ii) 2 +c I /’ 

1 0) J CO J 1 

arctan . = tan -1 {[(2+e ) 2 co 2 -4co 2 ]/JJn .oo 2 . } 

J u 7 j J 'j j 

- tan 1 { [ ( 2— c ) 2 co 2 -*lco 2 ]/^n . co 2 } . 

w J J J 

Consider now the interpretation of Eq . (7)- represents the 

band-limited time averaged power which is being accepted by the 

acoustic field inside the enclosure, that is, the power received 

over frequencies lying in band Aco . The band has a width Aco = 

c co where co is the center frequency of the band. c is a con- 
0) co 

stant percentage which defines the bandwidth (c = 0.232 for 

0) 

one-third octaves). The quantity < P^q >//c (JJ w as th e same quantity 
aS SpbfCw), that is, it is the power spectral density of the aver- 
age exciting blocked pressure which acts on the transmitting 
structure. <p 2 -^> as avera K e mean square blocked pressure as 

would be measured in band Aco obtained with a microphone sampling 
the surface pressures on the (hypothetical) blocked structure. 

As stated the transmitting structure has an area A exposed to 
volume V. The average surface mass of the structure is m. The 
structure has a mass law transmission coefficient t mt = (2pc /com) 2 , 
where p and c q are the density and sound speed of air. The nor- 
malization of structure modes is taken such that, generally, the 
modal mass is 

M = f m(x)ij j r ( x )dx = , 

A 

if the structure is unrestrained. The normalization of the 
acoustic modes is such that 
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The coupling of the exterior field to the rth mode of the trans- 
mitting structure is given in terms of the joint acceptance, 
j£(w). It is assumed that j£(w) changes little across Aw, i.e., 
that by evaluating j£(w) at the center frequency of the band, a 
good approximation is obtained all across Aw. For this reason 
Aw must not be too wide, typically one-third octave or less. If 
the exterior field is reverberant, j 2 (w) takes the appropriate 
form for that field. The same is true for a progressive wave 
field. 


The coupling of a structure mode to an acoustic mode of the 
volume is given by the term 

f ' (n , r ) = l I 4> n (x)i|j r (x)dx 


In the above x represents a point on the surface of the transmit- 
ting structure. Note that f'(n,r) is non-dimensional (hence the 
prime as a reminder) and that f'(n,r) is always squared. f ,2 (n,r) 
is a positive number less than unity. 

The term in the braces { } determines the coupling in frequency 

of the acoustic and structure modes (n and r respectively). The 

power flow depends on the locations of the resonance frequencies 

w^ and w^, i.e., relative to each other and to the band Aw, and 

upon the acoustic and structural loss factors n and n . 

n r 

Consider now the right-hand side of Eq.(l). W J4 is the band- 

diss 

limited, time-averaged power that flows out of the interior 
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acoustic field to the inner wall of the enclosure. This power 
flow is given by Eqs. (6) and (8) of Section III of [20] which 
reduce to 


W 


diss 


V 

pc; 


n to 2 

n n ^ 2 ^ 

<p > 

w n s , t 


n 


( 8 ) 


where <p 2 > . is the interior space average mean square pressure 

II u ^ b 

in band Am attributable to the nth acoustic mode. The development 
of Eq . (8) is presented in the appendix of [20]. 


The solution of a sound transmission problem Is obtained by 

setting the nth term of Eq. (7) equal to the nth term of Eq. (8), 

solving for the individual modal pressures, <p 2 > , , one-by-one, 

n s 

and by adding the results according to Eq . (3)> to obtain the 

interior level. 


3.3 General Results for the Noise Reduction 


The noise reduction of an enclosure is defined by 

<p? > . 

NR = -10 log 10 - 

e s,t 


(9) 


where <p 2 > , is the average mean square pressure over that space 

S S j z 

which will be occupied by the enclosure, and <p?> . is the space 

1 s 5 u 

average mean square pressure over that same volume when enclosed. 
Following the instructions given above for the solution of a 
sound transmission problem and letting, <p 2 > , = <p 2 >/2, it is 

6 S j L D X 

found that for reverberant field excitation, i.e., with 


J ‘ r 


Tsy, 
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ii n u) z 
n n n 


] 


arctan 


n 


c -c 
n r 


In + 
r 


t 


2c (b -b ) -b (c -c )" 
v n r r - 


■4n co 


r r 


3 -c n 
n r 


arctan 


( 10 ) 


Eq. (10) Is used in the present study to predict the noise 
reduction of the cylinder in one-third octave bands at those 
frequencies where it is necessary to account for the individual 
modal response and transmission contributions. At high frequen- 
cies, the summation over acoustic modes can be limited to those 
in the band (neAco). Eventually at sufficiently high frequencies, 
the acoustic modal density will become so large that Eq. (10) 
will require an excessive amount of computer time. To avoid this 
problem, Eq. (5) can be substituted into Eq. (2a) and evaluated 
according to the technique discussed in Section IV of 
Reference [19]. The reader is again referred to the paper for 
details. At these high frequencies, the net in-flowing power in 
band Ato is given by the expression 
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(ID 


which corresponds to Eq . (30), Section IV of [19]- 

In this expression t d is the resonant transmission coefficient 

n 

t R = (^ n r 1,Ic o r 'rad r 'iad /p “ 2mAi ir) ( 4pc o /A ) > 


where n^ is the structural modal density (modes/rad/sec) near 


and n r is the sum of the average external radiation 
i.e., ^pad* an< ^ ^he avera S e structural loss factor, n r . 


frequency co, 
loss factor, 

over Aco. Calculation of these quantities will be discussed later. 


The power absorbed on the wall is as before 


W 


diss 


V 

p c; 


n co 2 
n n 

to 


<p 2 > , 

s , t 


n 


Since the response is resonant acoustic, <o lies in Aco and the 
above reduces to 


W 


diss 


V 
pc ‘ 


o 


neAto 


p to 2 

n n <p 2 > 

to n s,t 


Further, noting that - co for any n, and letting the group of 
acoustic modes have an average loss factor n^, this reduces to 
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V 


. j — CO f] 

diss pc 2 n 


<D > 


neAu 


n s , t 


V 


wt O < P< > . 


pc^ ~ n s ,t , 


( 12 ) 


which is Eq • 7> Section III of [20]. Finally, setting 


ac S 
o 


(13) 


where a is the band average absorption coefficient and S the 
absorbing surface area, Eq . (12) reduces to the familiar form 


W 


diss 


aS 

4pc 

o 




( 14 ) 


The high frequency form of Eq . (10) is obtained by equating the 

results given in Eqs. (11) and (14). 


< t > e > s .t _ SS + 

.(V n rad) t R + T f] A 


(15) 


( X R + T f) A 

5 



where is defined to be the field incidence transmission co- 
efficient for mass controlled panel transmission, and is given 

by 


T 


f 


32 

it m 2 


I 

r<Aa) 


( 16 ) 
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Letting t_ + t_ = t, which is denoted as the total transmission 
n I 

coefficient, gives the form 


<P 2 > +. 

e s , t 



s , t 


ctS + 

T +( ^ > 

1 t r 

A 


tA 


(17) 


The noise reduction of the cylinder of concern in this study is 
computed using Eqs. (10) and (17). The frequencies over which 
the use of Eq . (10) is limited will be considered in a subsequent 

section of this report. 


3 . 4 Cylinder Noise Reduction 

Eqs. (10) and (17) are now specialized for the case of the test 
article, i.e., the unpressurized, unstiffened cylinder. Consider 
first the acoustic and structure modes, which can be ordered 
according to the occurrence of their resonance frequencies. 

For n = 1 , o> can be taken as the acoustic mode with lowest 
resonance frequency; for n = 2, the mode of second lowest fre- 
quency, etc. w can be similarly ordered. Corresponding to each 
acoustic mode is a trio of indices (q,n,s) and to each structure 
mode a duo of indices (M, N). These are used in the definition 
of the characteristic functions corresponding to each resonance 
frequency or eigenvalue. For the case of cylinder acoustics, the 
modal functions obey the Helmholtz equation 


subject to the requirement that the particle and surface velocities 
be equal along the surface normal 

9 $ 

W 2 - lkBcf, n = 0 ’ 


- 21 - 



where $ is the inner wall surface admittance (i.e., the ratio of 
surface velocity to surface pressure). The solutions for the 
mode shapes, <j> (£) = <b (r,0,z), and eigenvalues, A = A , 

for uniform surface admittance B<<1, are approximately the 
following : 


4) 


qns 


= J (m r) cos £n(0-0’)] cos 


Q TT Z 

L 


i = A - iA 
qns qns qns 


qns 


A 


qns 


0) 


— = / m 2 + (qfr/L ) 2 
c ns M 

o 


co n , e 

qns qns _ k qns 

2c 2A V 

o qns 


I 


Re [ g] I cf> 2 (x)dx , ( IS 
Q ns 


where q represents the number of axial half-waves, n gives the 
number of circumferential waves, and s counts within the sequence 
n of zeros of J’(ma) = 0. The geometry is as shown in Figure 1. 


The mode normalization obeys the relation 


e 


n 


£ 

qns 




(£)dV 


(19) 


For the cylinder 
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r- 9 2 ? 
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q n 


where 



2 q,n = 0 

1 q ,n > 0 


( 20 ) 


- 22 - 



FIGURE 1. CYLINDER COORDINATES AND NOMENCLATURE 



The nomenclature here may be confusing so to review, consider 
the meaning of Eq. (19)- This states that the mode whose reson- 
ance frequency is u> has a normalization such that the integral 
over the volume of the square of the mode shape is equal to the 
volume divided by a factor called which depends on the trio 
(q,n,s), and takes the value e which is given by the results 

C[T"1 S 

in Eqs. (19) and (20). Do not confuse n as a counter within the 
sequence w with the n in the trio of indices which defines the 
number of circumferential waves. The range of values for q,n,s 
is as follows: 


q 0, 1, 2, ... 
n = 0 , 1 , 2 , . . . 
s = 1, 2, 3, 


When values for q, n, and s are selected, the resonance frequency 
of the acoustic mode may be found using the relation 


O) 


qns 


o 


qns 


(rad/sec ) 


The remaining factors which are required for the cylinder calcu- 
lations using Eq.10 are the following: (1) the modal mass M^, 

(2) the coupling factor f'(n,r), (3) the joint acceptance 
j^fwT, and (4) the acoustic and structural resonance frequencies 
and loss factors. 


The structure mode shapes for a freely supported (ideal) cylinder 
are : 


<J> r (x) 


MN 

Hz,e) 


sin 


Mttz 

L 


cos N e 
sin N 0 


( 21 ) 
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The corresponding resonance frequencies are computed with the 
following characteristic equation [22] 


ft 6 -(K 2 +kAK 2 ) fj“ + (K 1 +kAK 1 ) ft 2 

- (Ko+kAKo ) = 0 , (22) 

where 

P r (l-k 2 ) 

ft = ^ a 2 ou r ,, T , k = h 2 /12a 2 

E I IjI 

P c is the mass density of the material of which the cylinder is 
constructed; E is the elastic modulus and y Poisson's ratio. Also 

K 2 = 1+| ( 3-y ) (N 2 +A 2 ) + k(N 2 +A 2 ) 2 

Kj= | ( 1-y ) £ ( 3 + 2y ) A 2 + N 2 + (N 2 + A 2 ) 2 

+ (n 2 + a 2 ) 2 ] 

K q “ | (1-y) [(l-y 2 )A 4 + k (N 2 + A 2 ) 4 ] 

with 

A = Mma/L . 

Constants K , Ki , and K 2 are from the Donnell-Mushtari shell / 
o 

theory; AK q , AKi , AK 2 are modifying constants for the theory 
of Arnold and Warburton [22]. 

AK 2 = 2 ( 1-y ) A 2 + N 2 

AKi = 2 ( 1-y ) A 2 + N 2 + 2(l-y) A 4 

- ( 2-y ) A 2 N 2 - | (3 + y) N 4 
AK 0 = |(l-y) [Ml-y 2 ) A 4 + 4A 2 N 2 + N 4 

- 2 ( 2-y ) ( 2 + y ) A 4 N 2 - 8A 2 N 4 -2N 6 J 
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Only resonances associated with flexure of the cylinder wall 
are retained. Using this cylinder model, it is found that 


m(2iraL)e 


M. 


where 


'N 


2 

1 


M 

r 

N = 0 
N > 0 


MN 


N_ mA _ 

r g n 


( 23 ) 


The interior coupling factor is 


f T (n,r) -*■ f'(qns,MN) 


•i-f 

A ’ 


, - * MN_ _ 
^qns x ^ ' x ' dx 


where A’ is the transmitting area. 


Two cases are considered, and are shown in Figure 2. In one case 
the cylinder and cavity lengths are identical; in the other the 
cylinder is taken shorter than the cavity. These cases were 
chosen after consideration of the application to the test article 
(the test cylinder configuration will be discussed in Section 4). 
For the case of Figure 2(a), 

L 2tt 


f ' (qns ,MN) = j 


n T / x , x qirz . Mttz Icos N0 

J (m a) cos n(0-0’) cos^ — sin-^ — { . ... } ad 

n ns L L Isin N0 


0dz 


o o 


f f 
qM 1 Ns 


0'=O 


f qM f 


Ns 


n = N, M q 
n = N, M = q 
n i- N 


(2b) 
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where 


J. T (m. T a) 
f _ N Ns 

Ns 2 N 


and 


f - -i- { 

qM 2tt ' 
= n 


1-cos ItflaliL i -cos C(M-q ) tt] l 

M+q M-q ) ’ q 


M = q 


Note here that 


f ' ( qns ,MN ) -► f ' (qNs ,MN) 

i.e., unless the circumferential waves are of the same order, 
the acoustic and structure modes are not coupled. This is true 
only for the case of an ideal cylinder, one in which there is 
uniform surface density. The actual test article has a lap joint 
running axially whose effect is considered in a later section. 
Note that it is possible to write Eq. (24) in the form 


f 1 (qns,MN) 


J N (nl Ns a) E N k I.l[ 1 - ( - 1>M< ' q ] 

2L ri_2 i_2 -i 


m q 
0 


n = N, and M / q 
n ^ N, or M = q 


where = Mtt/L, and = q-rr/L c . It can be seen that for Case (a) 
of Figure 2, f'(qns,MN) is zero if (M+q)/2 is an integer. 


Consider now the case of Figure 2(b). It is found that 


f ’ (qns ,MN ) 


f f 
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where 


f qM = COS 


COS TTC d (M-K) - COSH ( 1 + c )(M-k) 
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with 


k = k g L/n and c d = d/L . 


For the particular case where L = L, c = 0 and 

c d 


qM 2tt 


1 - CQSTT (M - k) ^ 1 - COSIT (M + k) 

(M - k) (M + k) 


f qM = o , 


which is the same as Eq. (24). 
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The joint acceptance function describes the coupling between 
the excitation field and the structure, and is defined by the 
relation 




A 2 S 


pbl 


(w) 


// 


Spbi (x | x ’ ,a))^ MN (x)if/ MN (x ' )dxdx' 


where A is the structure area, S ^(w) is the blocked pressure 
power spectral density, and S^^^(x|x',w) is the blocked pressure 
cross power spectral density. To evaluate the joint acceptance 
it is necessary to have representations for the excitation field 
and the structure mode shapes. 


Consider first the representation of the excitation field. It has 
been general practice, when representing random pressure fields 
with spatially decaying correlation, to use a correlation function 
of the form 

_ Re{S (x | x' ,(jo) } 

* sffl • 

p 

It is assumed in the present analysis that this function is 
separable in the longitudinal and transverse directions, i.e., 

p(I,u>) = p (£,a>)p (£,w) 
x y 


where £ = x'-x and c = y ’ -y . 


The representation for a diffuse excitation field is well estab- 
lished and is given by 


P x U,w) 

P y (?>w) 


sin (k^ ) 
k£ 

sin(k<; ) 
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For the case of the reverberant field, the joint acceptance 
functions for the cylinder must consider both sines and cosines 
of the circumferential wavenumber. One can write i- n the 

form 


•2 r< T\ 


2 rev 


,rev 


+ J MN Cu),4>) = j 




2^ eV x 

( CO 3 ijj ) 


rev 


+ J N 


>] 


where i p represents the cosine and <j> the sine of the wavenumber. 


The joint acceptance for the axial cylinder modes, 
given by Eq.(26) [23] : 
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Si and Cin are the sine and cosine integrals [24] 


The joint acceptances for the circumferential modes are [23], 
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where = y/2-na, y 2 = y'/2i\&, y and y' being circumferential 
coordinates of the two correlated points, and k the acoustic 
wavenumber. It follows that 


•2? e Y - 

4 ( m) = 


9 rev x 
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With this, all that remains to be determined are the acoustic 
and structure loss factors. These data are obtained experimen- 
tally and considered in Section 4. 


3.5 Transmission of a Discrete Tone 


Previously, it has been shown that for a noise field, the band- 
limited power flowing inwardly is 
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in 
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(29) 


This result is the same as Eq. (6), except j£(w) has been replaced 
by its defining relation. Note that only the cospectral density 
Cpbi (x | x ' ;w) , i.e., the real part of the cross spectral density of 
the exterior blocked pressure field, is used. 

For the present case, the external field is a propagating wave 
with a component along the z axis of the cylinder, with a discrete 
frequency co q . 
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Consider the cross power spectral density function 


c / — i — , v lim 
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Rp b ^ (x | x ’ ; t ) is the average cross correlation of the blocked 


pressure over the exterior, i.e.. 


R (x x 1 • x ) = — 
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-T/2 


Note here that the real part of p^^(x,t) is required. Consi- 
der the form of the blocked pressure field for the case of a 
cylinder with an incoming plane wave at angle <J> with 0 = it, as 
shown in Figure 3- TRe incident plane wave has amplitude and 
is given by 


’i = P i exp [" lu) o t + P ^ k x x + k z z) ] ’ 


where k = k sin <f>, k = k cos <)>, k being the acoustic wavenumber 
X z 

(a) o /c q . According to the geometry, x= r cos 0, which gives 
p^ = P^ exp ^-i(a 3 Q t - k^z ) + i k^ r cos 0 J 


This can be expanded in cylindrical harmonics using [25] 


00 

exp I + ik r cos 0| = /. e i“‘j (k r) cos mB , 

Lx J m m v x 


- 34 - 



z 



\ 


FIGURE 3. TONE TRANSMISSION: 
PLANE WAVE AT 9 = n 
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INCIDENT 




where 


1 m = 0 

,2 m > 1 


This gives the incident pressure wave as 
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The reflected (or scattered) wave can be written in the form 
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where H^^(k r) is the Hankel function of the first kind of 
m x 

integer order m, which represents an outward traveling wave for 

the assumed e la) o dependence. A are unknown coefficients which 

m 

remain to be determined. At the cylinder surface, the blocked 
pressure is 
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where p^ satisfies the boundary condition associated with a 
rigid (blocked) surface, i.e.. 
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which is equivalent to 
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Thus 


3 Pl \ 3p s 


3r A 3r 4 


Substituting and canceling out the t and z dependence gives 


P . e i m J ' ( k a ) 

1 z-t m m x 


cos me 


oo 

■-I 


A H ^ ^ ^ (ka) cos me 
mm x 


On a term-by-term basis, this reduces to 
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The blocked pressure Is therefore 
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At the surface x = (a,0,z). 
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Consider the term in the braces { }, i.e., 
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Now 
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Here use has. been made of an identity given in Morse and Ingard 
[26]. The term in braces becomes 
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Again 
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At the surface, the blocked pressure is therefore 
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Consider the limiting forms of the Hankel function and its 

derivative for k a>>m+l/2 
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Substituting this last result into Eq.(32) gives 
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Let 
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irk a 
x 


exp ^-i(u) Q t - k z z)J 


V cos m 0 -iY t -m-v 
x > — ^ e 'm(i ) 


m =0 


m 


But 


.m i 7 rm /2 

l = e j 
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so for the case 
P b l(x,t) 


+ TTm/2 


This result is consistent with Eq. (3.1.6) of [26], which is just 
a check of the results to this point. 


m=0 


cos m0 +i 

"E 6 

m 




m 


k x a>>m+2- 


4P, 


nk a 
x 


exp 


h 


i ( to t 
o 


V>] 


Consider now the average cross correlation of the blocked 
pressure over the surface 

T/2 

Rp bl (S|5';T) = Re[p bl (x.t)]Re 

-T/2 

The general result for the blocked pressure required here is given 
by Eq.(32). 


|p bl (x' , t+x )Jdt 


P bl (x,t) 


2P . 
1 

7T a 



k z z) ] 


x y e ^m+l cQs a ) 

m m x 

m=0 


The lagged pressure is: 


2 p . 

p bi ( x ' »t + T) = z^T 1 - ex P 


tt a 


[-U 


to t + 0) T- k z 
o o z 


■>] 


n=0 


e i n+1 cos n0'/H (1) '(k a) 
n n x 
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Since p bl (x,t) is periodic of the form e~ lW o t , the time average 
can be performed over one period, T q = 2tt/d o , giving 


R Pbl (x l x ' ; 


;T> = T n / 


V 2 


-T /2 
o 


Re [ p bl^ x,t ^J Re £ p bl ^ x ' >t+T dt 


Consider now the following theorem. If 

a = b = Se- la,t ; T = — , 


where A and B are complex numbers, then 


t f +T 

/ Re [a] Re ftl dt = | Re[A*B*] 

+- t 


P bl (x,t) = P bl (x)e _lw o t 


p bi(x',t + T) = P bl (x',T) e la) o t 


Rp bl (x |x , ;T) = | Re [P bl (x) P b *(x’,x)] 


II p 2 

1 D i io) t ik ( z-z r ) 

: x Re ' e o e z 

2 ( n 2 a 2 

00 a n+1 

^ e rn e n i cos mQ cos n ® ' 

" n 1 ^ ’ (k a) ’ (k a)* 

m,n=0 m x n x 
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This can be rewritten in the form 


Rp bl (x|x' ; t ) = ^ Re 


^ P i _i ^ 0 T+ k z ( z-z ’ 


2 2 
v a 


o 


\1 

[_m-0 


e cos m0 cos m0 ’ 
m 


H 


(1) T 


m 


(k x a) i z 


^ e^e^i (-i) n+ ^ cos m0 cos n0 ’ 

-s' n H^^’ck a) H (1) ’(k a)* 

m/n=0 m x n x 


The last term involving the summation of terms for which m^n is 
complex. Now note that the cospectral density, Cp bb (x |x' ;w) , 
required in Eq.(29) is to be integrated upon multiplication by the 
product ip (x)ip (x'), where 


,r/-x ,MN . Mttz , tq 

ip (x) = 4>(x) = sin — — cos N0 


\p r (x')= t|T(£t) = sln ^ 


MN 


Mttz ' 


cos N0 


Terms of the form cos m0 cos n0 1 will appear in Cp bl (x | x ' ;w) 
after Rp bl (x | x ' ; T ) has been transformed to obtain it. Since 

r 27T { 0;m / N 

I cos N0 cos m0 d0 = < tt;iti = N 

o ( 2 tt; m = N = 0 , 


it follows that all terms for which m / n vanish when the integral 

r,-v r / - , 


If 


X X 


Cp b i (x | x T ;oa) ip (x) p (x' ) dxdx’ 
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is evaluated. Thus, only terms for which m = n need be of 
concern and Rp^^(x|x';T) can be taken in the reduced form 


R Pbl ( x I x ’ ’ T ) = 


2P 2 

coslco T+ k (z- 

u 2 a 2 ^ ° z 


[“o T+ k z (z - z,) ] 


*2 


c e 2 cos m0 cos m0 ’ 


m 


m=0 


| H ( 1 ) ’ (k a) | 2 
m x 1 


Transforming gives the two sided cross spectrum 


Spbi ( x | x ' ; w ) = 


_ 1 _ 

2ir 


r 


Rp bl (x | x T ;t) e dx 


2 TF 


oo 

(&)I 

' / m=( 


c 2 cos m0 cos m0 ’ 
m 

*0 IhCD ’ (k a) | 2 

1 m x 1 


oo 

/ cos |co t + k ( z-z ' ) I 

L o z J 

m 


Now 


1 COT . 

e dr 


/ cos 

r CO T + k ( z-z ' )*] 

/ 

L 0 z J 


icoT . 
e dx 


0° 

j J e 1[nJ o T + k z (z - z ' )] + e - 1[ "o T + k z (z - z ' )] j e i U T dT 

i fe li y z - z,) 2 tt5(<i)+w 0 ) + e -ilc z^ s "" s 1 2ir6(w-o o )j 

= tt [,«.♦.„> e lk z (z - z,) + 6 Cto-w Q ) e- lk z (z - z,) ] 
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This gives 


Sp bl (x|x';u) = -r— \ 

TT^a^ 


P? e z cos m 0 cos m 0 T 

m 


m=0 


|h(D ' (k a) | 
1 m x 1 


x |^6(w+co o ) e^z^ 2 Z ^ + 6 (co-u) q ' 


The one-sided cross spectrum is defined by 


S pbl (x | x 1 ;cj) = 2 S DH1 (x | x » ; to ) ; 0 ^ u < 


Pbl 


2 P? e 2 cos m 9 cos m 9 ’ 


m 


2 2 
7T a 


m~0 


I H ( 1 ) ’ ( k a ) I 
1 ni x 1 


X e -lk z (z-z') 6 (ii)-w 0 ) 


Note that 

S Pbl^l^' ;w ' ) = S Pbl ^ f x ; uj ) 
Therefore, the integral 

Sp b i(x|x»;w) i|j r (x)ijj r (x ' ) dxdx' 

X X ’ 

is also equal to 

Spbl (x | x 1 ; to ) (x)^ r (x’ )dxdx’ 

X 1 X 




-ik ( z-z ’ ) 
e z 
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Letting 


^Pbl ( x I x > ^ ) Cp bb (x|x T 5 w) — iQp^-j-^Cx |x T jco) , 
it follows that 


II 


X X 


Qp^l (x | x ’ ; to ) ( x ) i^ r ( x ’ )dxdx; ' = 0 


i.e., the imaginary part of s p bl (x | x ' ;w) vanishes from the 
solution (Q Pbl = quadrature). Therefore, the only part of the 
cross spectrum to be used in Eq . (29) is the cospectrum 


2 P 2 " £ 2 cos m0 cos mQ ' 

Cp, , (x | x ' ;oj) = ) —7 cos k (z-z') 6(w-u) ) . 

Pbl ^a 2 L |H(D (ka)| 2 

m-0 1 m x 1 


Substituting into Eq.(29) and integrating with respect to w gives 


W 


2 tt /mA\ 2 x° 2_ 


in p v iTj ML tt o 


— «, e r\ 

V n n 

A 

n 

n 


Z frf J ' r) {/ / C pbl (x|x')<|/ r (x)ip r (x')dxdx'| 
n r r 1 x x f > 


GO ' 


"7 0) 2 V 2 T f/ 0)2 \ 2 T 5 

+ n n] l l ~ 4 ) + n r 


(33) 


where 

C pb i(x|x' ) 


2P? 

1 

•n 2 a 2 


m=0 1 m 


£ cos m0 cos m0 * 

— -7 cos k (z-z') 

|H(D (k°a ) | 2 z 
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with 


k° = k° cos d> 
z y 

T O ■, O • U 

k = k sin <p 
x 


o /2pc \ 

(“o /0 o) 003 * i T ML = 

(t» Jq/c 0 ) Sin <p 


The superscript o on k and k is to remind that the wavenumbers 

X z 

are evaluated at the excitation frequency w . Performing the 
integration, and noting that with r = (M,N), all cases m#N appear- 
ing in the summation m = 0 to 00 vanish, it is found that 



x x T 


^ r (x) 4> r (x T )dxdx ' 


L 2tt L 2tt„_2 2 i,„ „ n , 

2Pf Evt cos N0 cos N0 ’ 

IN . 0 , . v 

cos k ( z-z T ) 


fill 

J o J o 


iT 2 a 2 


|HU>'( k °a)p 


(3*0 


M 7T *7 M TT 7 ^ « 

x sin — — cos N0 sin — = cos N0 ’ a 2 d0d0 T dzdz' 


L 


L 


L L 


2P: 


4iT 2 a 2 


2 2 
tt a 


l H ] ^ 1) ' (k°a) ! 2 ^o ^0 


If 

J n J n 


, o , . v . Muz . Mttz ' , 

cos k ( z-z ’) sm — = — sin — = dzdz 

Z Li 


8P 2 L 2 


i H N 1),<k S a) i 2 


1 1 

// 

o o 


cos k°(z-z') sin Mttz sin Mttz' dzdz' 
z 


where 

k° = k°L ; z = z/L ; z' = z'/L 
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The integral is just the joint acceptance for a non-decaying 
progressive wave field, i.e.. 


-1 r 1 


J *(w ) = I I cos k°(z-z') sin Mirz sin Muz'didz' 
! '* Jr> Jn ' J 


In general, a progressive wave field with decaying correlation 
has the form 


j r > c ) 


f l f l r k zi i-z ’i" 

=j j exp — i — < 

J r\ l- rr -I 


, O | - 

k | z-z 

exp I Icos k° ( z-z ' ) sinMirzsinMirz * dzdz 1 

'o •'o L a z z 


For a = °°, this reduces to the non-decaying progressive wave 
z 

field. This is denoted by j^( ,a z =°°) . Eq.(33) therefore becomes 


W. 

in 


32P ! L2 /mAV o 5 V W , /T 

— (t) T ML% l 1%- ) n nJ 


* I 


f ' 2 (n,r ) J M^o ,a z °°^ 


r=MN 


M 2 a)“ 
r r 


id 1 ) ’ (k°a) | 2 
N x 1 


N) 2 + -*r 


(35) 


This is the fundamental result for the inflowing power. 


The power absorbed by the wall is 

, 2 


V 


1 i - i n n U) n 2 

abs pc 2 A w p n s,t * 


o 


n 


o 


where <p 2 > , is the space averaged mean square pressure attribu- 

n s 5 1) 

table to the nth acoustic mode. 
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Equating W in to W abg on a mode-by-mode basis gives (upon setting 
<p|> t = P?/2 = mean square pressure of incoming wave) 



In the above, as before 

e = £ .. = V/ f f [<p 2 , T dv 

n qNs ) ) J y qNs 

V 


///. 


V 


qNs 


(C)dv = 


m N a 2 -N 2 

■V J N (m Ns a ) 


2m 


Ns 


ttL 

2 


e e.t 
q N 


£ 


q’ 


'N 


2 ; q , N = 0 
1 ; q , N > 0 


and 


f ’ (n,r ) = f » (qNs ,MN) 


Also 


H N 1),(k ° a) - k x H N-1 < k x a) - ? 
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The range on N is 
M is 
q is 
s is 


0 , 1 , 2 , 

12 2 

0 , 1 , 2 , 

1 , 2 , 3 , ... 


Note, for instance 


h(D’ 

o 


(k°a) 

x 


k° H ( 1 ) (k°a ) 
x _! x 


- kVV<l>( k 0 a) 

X l x 


so 


l H o 1> ’ (k°a) I 2 


k° 2 1 J, (k°a) 
x 1 1 x 


+ i 


Y-, (k°a) | 2 


Also 


H n (1), (k°a) = k°H^ 1) (k°a) - 
1 x x o x 


i H a)(k°a), 


etc . 

Finally, the space average interior level at frequency to is 
found using 


<P 


2 

i 



V <P 2 > • 

Zj n s,t 


( 37 ) 
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The joint acceptance function required for progressive wave ex- 
citation is obtained from Wilby [ 27 ], 


with 






1 Oy 

. k L 

P M j l-(-l) e a cos k L 




T 0 T 

k L 
z 


o T , Mtt 


, ). / t \l f l _ “ . , U T . 1 V 1M * 

+ M-l) a z sin k z L + -y 



3 . 6 Response to Point Input Harmonic Mechanical Excitation 

The displacement response of the cylinder is 

W(x) = ^ 
r 


( 38 ) 
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£ r = ! r (w) is the modal displacement which is a solution of the 
equation 


-Y M 
r r 


+ pco 2 (I rr + J 


r rr 


W £ <I rs + J rs )5 s 
s^r 



(39) 


Eq. (39) is the same as Eq.(4), Section I of [ 19 ]. The above can 
be written in the form 


a £ 

rr r 


2 

r^s 


a 5 = 

rs^s 


-r 


pbl 


3 


or alternately as 


with solution 


Equivalently 


^• a rs-^s } { r pbl } » 


{ *r } = C5 rt ]{-r^ bl } 


^r 2 a rt r Pbl^ 

t 


SO 


W(x ,w) = -V i^ r (x) a r>t r^ bl (u) 

r,t 


(40) 
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t 

In the above Fp^^(co) is the generalized blocked force 


r Pbl^ = J P bl ( x (x)dx 

x 

For the present purpose, structural intermodal coupling will be 
neglected. This reduces the solutions for the modal displacements 
to 


l = -rn, , (w)/a 
r Pbl rr 


= -a r" 
rr P 


(a>) 


bl 


5 


giving 


W(x,w) 


- ^ ip r (x) a rr (a>) rp bl (w) 


(4l) 


where 


a 


rr 



M 


r r 


+ 


2 

pu> 



Now consider the generalized blocked force. For the case of point 

harmonic mechanical excitation at x = x = (a,0 ,z ). 

o 5 o 5 o 5 


P bl (x,w) = F(oj) 6(x-x q ) 

Note that p b ^(x,t) = F6(x-x Q )e , i.e., F(u)) is the amplitude 

of the force when the frequency is w. 
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In this case, for the cylinder 


6(x-x ) = — 6(0-0 ) 6(z-z ) 
o a o o 


Therefore (for the ideal cylinder) 


r 


r 

Pbl 


(w) 


nL 

F(gj) 

a 


6 ( 0-0 q ) 6 (z-z q ) 


x sin cos N(0 -0j)ad0dz 

Mirz 

= F(u)) sin —= — — cos N(0 -0, ) . (42) 

L o 1 

0-^is a reference angle which is set equal to 0 Q for the case of 
the ideal homogeneous cylinder. This gives 

Mtt z 

r r / n v , o r MN , s 

r pbl (uj) = F ( (JO ) sin ^j- — = r p bl ( w ) 

Substituting, the response is found to be 

W(x,u)) = - sin cos N(0-0 q ) 

M,N 

* { -V»N + <»* (l™™ ♦ 

Mttz 

x F(co) sin — z - — ■ 
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The acoustic effects are negligible in the sense that the acoustic 
reactance is much smaller than the shell reactance. The internal 
radiation resistance is small compared to the external which may 
or may not be comparable to the mechanical resistance. Thus we 
can write 


Y M 
MN MN 


+ pa)' 


6 


MNMN 


jMNMN ^ 


as 


where 



in which is the sum of the mechanical (structural) and 


sin 


external radiation loss factors. 

Mttz^ 

oo • 

C7 “I V'l 

W(z,0,w) = F(oo) 

M-l 
N=0 


Therefore , 
Mttz 


Y M 
MN MN 


sin 


L cos N(0-0 q ) 


(^3) 


The time mean square displacement is 


<w 2 (z,0 )> t 


W(z,Q,aQ 1 2 _ 
2 



Therefore , 
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<w 2 (z , 0 ) > 


= j Re 


1 1 

HN mn 


Mttz niTZ 

sin — ^ — — sin —= — — 


Y Y M M 
MN mn U MN mn 


x sin cos N(0-0 ) sin — cos n(0-0 ) 

L o L o 


] 


Mttz 

2 ^ sin 2 — T — — 

sin 2 H12L5. cos 2 N( 0-0 ) 

Y I 2 M 2 L ° 

• l mmI 


2 


MN 1 MN 1 MN 


• r».z r 


Mttz ttittz 

o . o 

in — = — sin 


j_, L . Mttz . . , n n \ 

y y ^ 5w? jyr SIH T” COS N(B 6 _ ) 


MN mn MN mn MN mn 
MN i mn 


miTZ r n „ \ 

x sin — — cos n(0-0 ) 
Lj • o 


Consider the summation in which MN ¥■ mn. For every pair of 
modes r = MN, s = mn, there is a corresponding pair r = mn, s = MN 
The sum of the two is 


Mttz 


muz 


sm- 


-sm- 


o . Mttz . ttittz Q a \ 

sin— r — sin— — cosN (0-0 Q )cosn(0-0 o ) 


L 


M,.„ M 


mn 


1 


1 


Y y * Y Y * 

MN mn mn MN 


This can be written as 


Q.... ( 0 , z ; 0 , z ) 

MNmn ’ ’ o o 

M M 
MN mn 


W~ T 5? 

Y Y Y Y 

MN mn mn MN 


20.... ( 0 , z : 0 , z ) 

MNmn * 5 o’ o' 

M M 
1 MN mn 


RetY MN ]Re[ V ] + Ira[Y MN ] 


MN 1 


mn 
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which is a real quantity. Thus it is possible to write the time 
mean square displacement at a given point on the cylinder in the 
form 


<w 2 (z,0 )> t 



Q MNMN^ 6 ,Z;e o >Z o :) 
I y mn 1 2 m mn 


+ 


2 I 

MN mn 
MN / mn 


0 „ ( 0 , z ; 0 z ) 

V MNmn ’ o o 

Y Y* M M 
MN mn MN mn 


(4H) 


v.'here 


0 ( 0 , z ; 0 , z ) = sin — r 

w MNmn o’ o' L 


Mttz mTrz 

o o Mttz . muz 

sin — sm — ^ — sm — = — 


L 


. cos N(0 - 0 ) cos n(0 - 0 ) . 

o o 


This is a fundamental result. 
The acceleration is 


<a 2 (z s 0)> t = (jj'*<W 2 (z,0 )> . 

The space-average mean-square displacement is 


<w 2 (z,0)> t - 


1 

2iraL 



< w 2 (z,0 )> t 


ad0dz 
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Since all the cross terms integrate to zero. 


<w 2( z,0 )> t; j 


1 f f F 2 V Q MNMN^ e,Z;0 o’ Z o 

2 W 0 l 2 ^ 1 y mn I ZM MN 


F 2 

¥ttL 


2 


sm 
f MN 


Mttz 

2 O 


L 

"Y 2M2 
1 MM 1 


MN 



sin 


2Mirz 

L 


dz 


x 



cos 2 N( 0-0 o )d0 . 


Now 


where 



'N 


2 N = 0 
,1 N > 0 


Therefore , 


<w 2 ( z,e)> t; - 


F 2 

T 


I 


C N Sin 

I y mn 1 2 


Mttz 

o 
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ad6dz 



For the cylinder 


M 

MN 


£ N mA _ £ N M T 

4 


where M T is the total mass of the cylinder wall, so 

Mttz 
2 0 

L 

y r 5 " * 

MN 1 


<w 2( z,0)> tj - = 


MN 


sin 


'N 1 


( 45 ) 


This is a fundamental result. 

Interior Sound Level - Point Harmonic Mechanical Excitation 
At any point £ in the cavity, the pressure is 


p 2 (C) = -pm 2 / G (x | £;co)w(x)dx 


i 


( 46 ) 


where the interior space Green's function is 


4> (xH (O 

G (x|£;w) = Y - — 

P n ///4> n (C)dv(A n ~ k ' } 


Substituting the results for the cylinder gives 


■ 2 

q,n,s 


G p (x | £;w) 


c cos n(0 
qns 


G (a , 0 , z | r • ,0' , z * ;w) 


0')J (m a)J (m r')cos — cos 
n ns n ns 


qnz 

L 


qiTZ ' 
L 


V(A 2 - k 2 ) 
qns 


( 47 ) 
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According to Eq . (43), for harmonic force F applied at z o ,0 q 


W( x ) = W ( z , 0 ) = F 


1 


M ttz q 

sin T, *** • Mttz tm/a a \ / it o \ 

y — jvj sin -y— cos N(e~e ) - (^8) 

Tl/TKT^ /T > T -Ll ^ 


MN MN MN 


Substituting Eqs.(47) and (48) into Eq. (46) gives the interior 
pressure as 


Mttz 


e F sin — j- 

qns L 


P n -(I) = Pj (r * , 0 ' ,z ' ) = -pm 2 Y A 

^ m v(I 2 -k 2 ) y mm m 


o 


qns ii. 


qns 


MN MN 


x J (m a) J (m r') cos 

n ns n ns L 


2tt 

x | cos n(O-O') cos N(0-0 ) ad0 

to 


l 

f L < 

x / sin 

Jo 


Mttz 


qtTz 

L 


Now consider the trigonometric identity 

cosn(0-0') = cosn0 cosnO' + sinn0 

The product required under the integral 



sinn0 

above 


is 


cosn(O-O') cosN(0-0 Q ) = 

cosn0 cosNO cosnO' cosN0 + cosN0 cosN0 

o o 

x sinn0 sinn0' + cosnO cosnO' sinNO sinN0 o 
+ sinnO sinn0 ' sinNO sinNO 

o 
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In the integral, all cases in which n / N vanish. 


Also , 


1 


2 TT 


COS N0d0 = TT £ 


N 


I 

l 


sin 2 N9d0 = 


tt N > 0 

0 N = 0 


cos N0 sinN0d8 = 0. 


Therefore the integral 

-2 TT 


f 


cos N( 0-0 T ) cos N(0-0 Q )ad0 


= e^iracos N(0 T -0 o ) 


Nov/ let 


D. 


Mq 


1 


L 


. Muz qirz , 

sin — — cos dz = 


Lf qM M * « 
0 M = q 


Then the Interior pressure becomes 


Mttz 


P 1 (r’ ,9’ ,z* ) = 


2 f 

■pm 7 


'qNs 


sm 


o 


i—l ( T2 _y.2 \ 

MNqs ' A qNs K ; 


M Y 
IN MN 


J N (m Ns a ' )D Mq 


qw z 1 


x e j^TTa cosN( 0'-0 q ) J N (m Ns r') cos L 


(^9) 


This is a fundamental result. 
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The space-averaged mean square pressure is 


<p2 (r T > 0 ' , z ’ )> Sjt 


-§y P* ± (r ' ,0 ' ,z ' )] r ' dr ’ 

v 


* d0 ' dz ' 


Consider the term 


f/J[ COS M(e '" e o ) W’’ cos nr] 


V 


x jjcos N ’ ( 0 1 — © Q ) 'Tjj, ( m N* s » r ' ) cos ^ T p Z * j r ’ 1 d6 ’ dz * 


This is non-zero only when N=N' 5 s = s', q=q ’ and becomes simply 


Iff 

V 


Vs (?)d ' ; = V/e qNs 


The space-averaged mean square pressure in the interior therefore 
reduces to 




p 2 c 4 a> 4 A 2 

2 F 2 

2V 2 


V 4 Mirz 

x /, sln — r 2 " 2 (n,r)-i 

r=M,N 


L qNs 


n=q,Nj s 


3 * / o 

M 2 


CO 


n 


[ (1 -^ ) + , 1 n] “r[ (1 -^> + 0 r] 


(50) 
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3 . 7 Interior Response for Progressive Wave Field Excitation 


The solution for the band-limited case can be obtained by equating 
W^ n of Eq . ( 7 ) to W^ ss of Eq . (8). Here, as opposed to the case 

for the noise reduction where the ratio of the interior pressure 
to the exterior diffuse level was sought, the ratio of the space- 
average interior to the average blocked pressure is to be deter- 
mined. Eqs . ( 1 ) , ( 3 ) j (7 ) , and (8) give the solution to be used 
at low frequencies. 


<P?> 


i s , t _ 4 A 2 / mA 


<p bl > s , t 


C 03 T T 2 
03 V 




I •» 2 


n 


j 2 ((o)f' 2 (n,r) 

M 2 D 
r nr 


x 


+ 



2c (b -b )-b (c -c ) 
n r n n r n 



arctan 

n 

arctan 1( 51 ) 


The result is similar to Eq . (10), except that now j£(w) (=j^( w) ) 

is taken as the joint acceptance for the progressive wave field. 


The joint acceptance is considered to be separable in the axial 
and circumferential directions, i.e.. 


J mn ^ u ^ J ]v[ ^ w 3 j ^ ( to ) 
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The axial component, j^(<o), is given by Eq . (26), where a z is 

allowed to take on finite values to account for the decay in 
correlation. The joint acceptance for the circumferential bend- 
ing modes is [23] 



where k is the circumferential trace wavenumber, L = 2na is the 

circumference of the cylinder, and a^ is the decay coefficient. 

The decay constants a and a , and the trace wavenumbers k and 

y z y 

k usually must be measured for any excitation field configura- 
tion . 


The high frequency result is obtained by equating the net inflow- 
ing power given by Eq . (29) of Reference [19]j i.e., Eq . (53) 
below. 
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(53) 


to the power dissipated given by Eq. (14). The internal radia- 
tion resistance required is found with the relation 
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(54) 


T3 int ( \ 

R rad (M) 


2pc 


k 2 A 2 


o 


rev 


where < > r indicates an average over all structural modes 
resonant in the band. 

3 . 8 Modifications for the Non-ideal Cylinder 

Figure 4 shows the unpressurized, unstiffened cylinder that was 
used as the Phase I test article. The cylinder differs from an 
ideal cylinder in that there is a lap joint running longitudinally 
0.00457 meters ( 0 . 1 8 inches) wide. At the joint, the cylinder 
thickness is doubled, as shown in Figure 5» This gives the 
cylinder a non-uniform surface density. Because of the narrow 
width of the joint in comparison to the cylinder circumference, 
the joint can be modeled as a mass discontinuity located at 
arbitrary angle 0=0. Previously it has been noted that in 
the case of an ideal cylinder, acoustic modes with circumferential 
orders n can be excited only by structure modes with circumfer- 
ential order N = n. The lap joint complicates this simple coup- 
ling, causing all of the acoustic modes to be excited by all of 
the structure modes, as will be shown. 

In addition to the above, the test article has end caps which are 
not rigid. While the noise reduction capability is high, the 
caps do flank at extremely low frequencies, i.e., below the first 
acoustic mode of the cavity. 
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FIGURE 4. PHASE I TEST ARTICLE: U N P R E S S U R I Z E D , UNSTIFFENED CYLINDER 





FIGURE 5. 


CYLINDER LAP JOINT DETAIL 
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3.8.1 Ci^Vinder Lag Joirvt Ef_f_ect_ 


In general, the cylinder wall response for an arbitrary harmoni- 
cally applied load is given by 


w(x) 



;o> )q ( x T ) dx ' 


(55) 


where q(x') is the pressure load, G(x|x T ;w) is the cylinder 
Green's function. 


and 


G(x | x* ;oj) = 


I 


(x)ip r (x') 

1T7 


r r 





Consider the case where q(x r ) is a line load along the cylinder 

at 0 = 0 , i . e . , 

o 5 5 


q ( x ' ) 


F^ 

ci 


6 ( 0 '- 0 o ) 


q ( z ' ,0 ' ) 


(56) 


F* is the force per unit of length 


in the z direction (axially). 


Let the damping be zero, i.e., take n r = 0. Then Y^w) = . 

The question to be answered is how do the mode shapes and resonance 
frequencies change if F' is due to the vibration of a mass dis- 
tributed along ( z , 0 ) . 

In terms of the cylindrical coordinates, Eq.(55) becomes 

n L 

G(z , 0 | z ' , 0 ' ;m) q(z ' ,0 ' )ad0 'dz ' 
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Now let 

P' = -w 2 m ’ ( z ’ ) w 


where m* is the mass per unit of length of the line load, and w 
is an arbitrary amplitude of motion (meters). The hypothesis is 
that as the frequency of motion tu approaches to (1—6 ), w(z,0) 
takes on the value ’t' r (z,0) - ij; r (z,0) which is the new mode 
shape accounting for the mass discontinuity, provided that m* 
is not too large. Also, w^d-6^) is the new resonance frequency. 
The negative sign is chosen because of the mass being added. 

& r is a constant, presumably small for the case of the lap 
joint, which must yet be determined. Substituting for p* gives 
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In the above, lower case m and n are used to distinguish the 
mode mn from the mode MN. An approximation for the mode shape 
for the mnth mode considering the lap joint is 
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From this form. It is seen that one can take 


<f> ( z ) = \p ( z ) 
r m r m 


and write 


r m ' w ^ n ( 9 0 ) 


V 9) = 


r 2 it 

2m6 / if) 2 ( 9 )ad0 

mn / r n 

Jo 


ip (0) 
r n 


m’w V 

m 


03 2 li) » T ( 0 )lk T (0) 

mn r N o r N 




CO ^ "I /" 

1-^2. / 

^mN-Uo 


2 IT 


if)^(0 )ad0 


Now multiply through by 


m ' will (0 ) 

r n o 

2m6 f 2ir if) 2 
mn | r n 

Jn 


" 1-1 


(0)ad0 


- 72 - 


and let 



^(e)ade 


= y 


n 


Renormalize and call the new mode shape Y (0)(note the capital 
letter) by setting 
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(57) 


This result states that the lap joint causes the circumferential 
mode shape i|; ( 0 ) of the ideal cylinder to take the new value 
V (0); this mode shape is approximately equal to 4^(0), but 
differs in that components of all other n order are present. 

Eq. (57) may be thought of as a Fourier representation of ¥ n (0) 
over the range 0 ^ 0 ^ 2tt, i.e.. 


where 
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As can be seen, the values of the coefficients A^ depend on the 

shift in resonance frequency caused by the addition of the line 

load. Thus, to solve for the coefficients A^, an estimate for 

6 is needed as a function of m'(i.e., the added mass per unit 
mn 

of length). It can be shown that 
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where B is the overlap cross-sectional area. A, divided by the 
radius of the cylinder, i.e.. 


B = A/a. 


This gives 
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3.8.2 o/ End Cags^ and Lon£itu^J 1 iml FlexibiJAtj^ 

The test cylinder has a wall thickness of 0.0016 meters ( 0 . 0 6 3 
inches). The end caps are 0.0127 meters (0.050 inches) thick. 
Nevertheless the noise reduction of the test article can be 
limited at low frequencies, i.e., below the lowest acoustic mode 
of the cavity, by the end caps. The problem is illustrated in 
Figure 6. At low frequencies, the noise reduction afforded by 
the cylinder wall can be high, higher perhaps than the noise 
reduction afforded by the end caps. Whenever this occurs, the 
effects of the end caps must be taken into account. 

The cylinder also has longitudinal flexibility. Whenever a uni- 
form pressure acts on the cylinder, the cylinder shortens. This 
occurs in addition to the bowing-in of the end caps and shrink- 
ing of the circumference of the cylinder. The shrinking of the 
circumference is accounted for in the dynamics of the present 
cylinder model. However, the other two effects are not included. 
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Low Frequency - Long Wavelengths = Uniform Exterior Pressure 








Following Lyon [14], the exterior pressure p 0 produces an 
interior pressure p^ through a volume deformation of the end-cap 
denoted by AV. The interior pressure is p^ = k c AV where k is 
the stiffness of the acoustic volume, i.e., k = 1/C , C being 
the acoustic compliance = V/ pc 2 . The volume displacement AV 
is determined from the pressure difference acting on the end cap, 
q = p-p^, and the end cap compliance, C p , i.e., AV = (PgPpCp* 
But AV is the same by both calculations, hence 
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This gives 
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and the noise reduction as 


NR = 20 log (l + ^ ) . 


Hence an expression for the compliance of the end-cap is required 
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Roark [28], gives the displacement of a circular plate, of radius 
a, with simply supported edges, loaded under uniform pressure q, 
as 
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and 


G 11 = 1/64 

q is the load per unit of area in psi, and D is the bending 
rigidity, D = Et 3 /12 (1-p 2 ) . 


The volume displaced per unit pressure is 
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(59) 



Longitudinal Flexibility 


If the cylinder stiffness as a column is taken into account, it 
is found that the compliance, (V = AV/q is 


r » - rca 3 L 

p 2hE 


( 60 ) 


Test Article Noise Reduction 


The noise reduction owing to the end caps and longitudinal flexi- 
bility is denoted here by a bar, i.e., 

™ / C \ 

NR = 20 log 1 + 


C + C' 

P P 


( 61 ) 


where C p is given by Eq . (59) and by Eq.(60). 


If NR is the cylinder wall noise reduction, then the test article 
noise reduction, NR^ , is 


NR = -10 log 


1 ( 10 -NR/10 + 1q -NR/10 


) 


( 62 ) 
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4.0 EXPERIMENTAL INVESTIGATIONS 


The Phase I study has a two-fold objective: development and 

validation of preliminary analytical models. The development 
of the analytical models was presented in Section 3- Section 5 
will consider the validation. In order to validate the models, 
an experimental program had to be undertaken to obtain a reason- 
ably comprehensive data base. Experiments were done in joint 
effort at NASA Langley Research Center, under the supervision of 
Mr. William H. Mayes and Mr. Conrad Willis, and at BBN in Los 
Angeles. In some cases, identical tests were performed at both 
locations. The series of tests was the following: 

1) Measurements of the structural and acoustic loss 
factors of the cylinder and of the cavity. 

2) Measurement of the noise reduction of the cylinder, 

i.e., the difference between the exterior and the interior 
sound pressure levels given a random, diffuse (reverber- 
ant) exterior acoustic field. 

3) Measurement of the interior sound level given excitation 
by a harmonic plane acoustic wave incident at various 
angles relative to the cylinder axis. 

4) Measurements of the interior sound level and the 
acceleration response of the cylinder given a point 
harmonic (tonal) mechanical excitation of the cylinder 
wall. [The difference between the interior pressure level 
(or the cylinder acceleration level) and the applied force 
level was measured using a broad band noise force input as 
opposed to a tone. The resulting "transfer functions" are 
directly comparable to the tonal predictions.] 
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5) Measurement of the interior sound level for the case of 
excitation by a random, progressive wave field with 
decaying correlation. 

The tests related to the measurements of the acoustic and 
structural loss factors and the cylinder noise reduction were 
performed both at LaRC and at BBN . The tone transmission, 
mechanical excitation, and progressive wave field experiments 
were done at BBN. In addition to the primary tests above, 
measurements were made of the cylinder modal characteristics for 
comparison against the analytical model. 

4 . 1 Phase I Test Arti cl e 

The cylinder of concern has been briefly described in Section 3 
(refer to Figure 4). The cylinder is 1.2192 meters (48 inches) 
long and has a wall thickness of 0.0016 meters (O.O 63 inches). 
Inserted into each end of the cylinder is a 0.0127 meter (0.50 
inch) thick ring, which has an inner diameter of 0.4064 meters 
(16 inches). The rings shorten the cylinder somewhat relative 
to the length of the enclosed volume. They also disturb the 
otherwise perfect nature of the cylindrical cavity. However, 
the effect on the cavity shape is not considered important. 
Basically, the end rings affect only the cylinder structural 
dynamics, i.e., they make the cylinder shorter by stiffening 
the ends. At the same time, the rings provide the important 
function of maintaining the ends circular. The cylinder is 
considered to be 1.1938 meters (47 inches) long because of the 
rings . 
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A small acoustic source was mounted Inside the bane cylinder 
and driven with broad-band noise. A Bruel and Kjaer (B&K) Type 
4131 25mm microphone was used to measure the resulting sound 
level and the corresponding spectrum was displayed on a Spectral 
Dynamics SD 360 analyzer in order to identify the resonance fre- 
quencies of acoustic modes. Figure 7 shows a typical response 
spectrum, in which the lower-frequency resonances have been 
identified. Acoustic resonances up to about 850 Hz have been 
identified by comparing measured and calculated results. These 
are shown in Table 1. 


Figure 8 shows the spectrum of the acceleration response of the 
cylinder surface given acoustic excitation from inside the 
cylinder. Many peaks in the spectrum correspond to the acoustic 
modes of the contained volume. The others are the vibration 
modes of the cylinder structure itself. Some of these are 
identified in Figure 8 and a more comprehensive comparison is 
shown in Table 2. Fair agreement exists between measurements 
and the prediction models, as expected for this ideal configur- 
ation. It is clear that as frequency increases, both structure 
and acoustic modes accumulate quickly and are so closely spaced 
that their resolution becomes impractical — hence the recourse 
to statistical analytical models whenever appropriate. Note 
that the resonance peaks for both structure and cavity are very 
sharp, characteristic of a very lightly damped system. More 
heavily damped systems have wider modal bandwidths and resolu- 
tion becomes more difficult still. 
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FIGURE 8. CYLINDER ACCELERATION RESPONSE TO INTERIOR ACOUSTIC EXCITATION 





Table 1. Acoustic Resonance Frequencies of LaRC Cylinder 


Mode Order 
q N s 

Calculated Frequency 
(Hz) 

Measured Frequency 
(Hz) 

1 

0 

1 

140. 7 

142.6 

2 

0 

1 

281.4 

282.8 

0 

1 

1 

395.7 

400.0 

1 

1 

1 

420. 0 

— 

3 

0 

1 

422.1 

426.0 

2 

1 

1 

485.6 

490.8 

4 

0 

1 

562.8 

568.0 

3 

1 

1 

578.6 

585.0 

0 

2 

1 

656.4 

641.0 

1 

2 

1 

671.3 

670.0 

4 

1 

1 

688.0 

685.0 

5 

0 

1 

703.4 

709.0 


L c = 1.2192 m 

Cavity Dimensions a = 0.25** m 

c Q = 3**3 m/s 
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Table 2. Structural Resonance Frequencies of LaRC Cylinder 


Mode Order Calculated Frequency Measured Frequency 

M N (Hz) (Hz) 


1 

4 

12 6.6 


131.0 

1 

3 

153.4 


153.5 

1 

5 

159.7 


163.0 

1 

6 

221.9 


220.0 

2 

5 

264 . 0 


264.0 

2 

6 

272.2 


275.0 

1 

2 

289.0 


283.0 

1 

7 

301.1 


297.5 

2 

4 

331.0 


318.0 

2 

7 

327.8 


322.5 

3 

6 

400.7 


387.0 

3 

7 

403.0 


403.0 

1 

8 

39^.1 


407.0 

2 

8 

410.9 


424 . 0 

3 

5 

456.9 


453.0 

3 

8 

456.6 


464.0 

1 

9 

500.0 


479.0 

2 

3 

510.0 


506.0 

2 

9 

512.6 


529.0 

4 

7 

531.6 


539.0 


Cylinder 

Dimensions 

L = 1.1938m (47 in) 
a = 0.254m (10 in) 
h = 0. 00160m ( . 063 in) 
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Figures 9 and 10 show the computed resonance frequencies for the 
cylinder and cavity over the range of N orders required for the 
present study. The cylinder data are for the 1.1938 meter 
(47 inch) cylinder and differ only slightly from the data for a 
1.2192 meter (48 inch) cylinder. 

4 . 1.2 Damping 

Estimates of structural and acoustic loss factors are required 
for use as input data for the various analytical formulations. 

In general the loss factors vary from mode-to-mode . When band- 
averaged measurements are made, the resulting loss factor is 
determined by the most lightly-damped, most strongly excited 
modes. In general, since the modal forcing function for the 
driven modes can vary significantly over the band, different 
excitations can produce different values of band-averaged loss 
factor, which means that the actual excitation under investi- 
gation should be used in the measurement of the average loss 
factors . 

Measurement of the decay of steady-state vibration energy levels 
after removal of the excitation gives correct estimates of the 
average loss factors for mechanical or acoustic excitation, but 
when the coupling, for example, controlling the power flow from 
the vibrating cylinder to the interior acoustic field depends on, 
say, only one out of a group of modes resonant in the frequency 
band of interest, large errors can be produced. In such cases 
detailed modal measurements are required, and this can be done 
using single-frequency excitation tuned to the particular modal 
resonance frequency of interest, measuring the decay rate of the 
modal energy level or by using a broad band noise and measuring 
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FIGURE 9. CYLINDER STRUCTURE RESONANCE FREQUENCIES 
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the bandwidth of the resonance. Both approaches were employed, 
particularly at low frequencies where modes were widely separated 
in frequency. 

The "response-decay" method involves measurement of the rate of 
decay of vibration or pressure level from its initial steady- 
state level after the excitation has been removed. Figure 11 
shows the instrumentation used. The amplified response signal 
is passed to a linear-log converter and then to a variable 
time-base storage oscilloscope, where the decay rate (ideally, a 
straight line) is measured with a calibrated protractor. For 
"response-amplitude" measurements, the analog signal for broad- 
band excitation is passed to a Spectral Dynamics SD 36 O Narrow 
Band Analyzer, fitted with translators, and the bandwidths of 
particular modes measured by fitting a smooth curve through the 
digital output response data close to the resonance frequency. 

Acoustic loss factors were measured using both methods out- 
lined above, and close agreement was found for the resonant modes. 
Data were also taken using the response decay method at frequen- 
cies well removed from the resonance peaks, thereby reflecting 
the absense of resonant response. The loss factors were generally 
an order of magnitude higher than those measured on or close to 
the resonance frequencies themselves. The measured data are shown 
in Figure 12. It is clear that the loss factors of the resonant 
acoustic modes of the bare cylinder are very low. 

The LaRC data are also shown in Figure 12. The loss factors of 
concern correspond to the LaRC case measured without the lead 
cover. The LaRC data appear to give the damping characteristics 
of the nonresonant acoustic response at low frequencies. 
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FIGURE 11. 


TYPICAL INSTRUMENTATION FOR LOSS 
FACTOR MEASUREMENT 
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Frequency, Hz 


FIGURE 12. ACOUSTIC LOSS FACTORS OF BARE CYLINDER 
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A precise explanation for the measured behavior, i.e., low 
modal loss factors for resonant response of the modes, and high 
loss factors for nonresonant response of the modes cannot be 
given. It is believed to be caused by different loss mechanisms 
(absorption and radiation). The nonresonant acoustic response 
is composed of many acoustic modes that may be driven above or 
below their resonance frequencies. As a group, these modes have 
more ways to shed their energies; whereas a single acoustic mode 
excited at or near resonance has not so many options. However, 
because of the fact that the response is nonresonant, the dis- 
sipation should play little role in determining the response (the 
response will be of a forced nature, either stiffness or inertially 
controlled). Thus, the average loss factors for the one-third 
octave bands where no modes are resonant should not influence the 
interior predictions. Knowledge of the loss factors for those 
bands where acoustic modes will respond in a resonant fashion is, 
however, very necessary. Table 3 summarizes the data which were 
used as input to the computer program. Where available, the 
bandwidth measurements were averaged to obtain a value for each 
band. All modes resonant in that band use the average value 
during computations. Where no acoustic modes are resonant in a 
band, the acoustic loss factor corresponds to that measured with 
the response decay method for tonal excitation. 

Structural loss factors were measured using the response decay 
method with magnetic, acoustic, and point impact excitations. 

Modal data were also obtained using the bandwidth method 
described previously. The measured data are presented in 
Figure 13. Also shown are LaRC measurements. In general, the 
structural loss factors lie in the range below 5 x 10“^, which 
of course is indicative of a fairly lightly damped system. 
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Table 3. Band Average Loss Factors--Undamped LaRC Cylinder 


Frequency (Hz) 


Acoustic 
Loss Factors, q 

n 


Structural 
Loss Factors 


50 

0.143 

— 

63 

0.143 

— 

80 

0.120 

— 

100 

0.107 

— 

125 

0.0031 

0.0030 

160 

0.100 

0.0035 

200 

0.155 

0.0027 

250 

0.050 

0.0045 

315 

0.0087 

0.0023 

400 

0.0019 

0.0015 

500 

0.0014 

0. 0017 

630 

0.0044 

0.0025 

800 

0.0025 

0.0036 

1000 

0.0097 

0.0018 

1250 

0.0065 

0.0013 

1600 

0.0032 

0.0020 

2000 

0.0042 

0.0025 

2500 

0.0025 

0.0023 

3150 

0.0019 

0.0025 

4000 

0.0013 

0.0019 

5000 

0. 0011 

0.0020 


n 


r 
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FIGURE 13 


MEASURED STRUCTURAL LOSS FACTORS FOR 
BARE CYLINDER 







Table 3 gives the structural loss factors used for computations. 
If a structure mode is resonant in a given band, it is assigned h 
loss factor equal to the value listed in Table 3 for that band, i 
the same manner as is done for the acoustic modes. The structure 
data in Table 3 are a composite of the BBN bandwidth and impact 
measurements. There are no structure modes resonant below 125 Hz 

4 . 2 Noise Reduction Measurements 

Noise reduction measurements were conducted at NASA LaRC and at 
BBN on the bare cylinder. The excitation was a random, diffuse 
(reverberant) acoustic field. The layout and instrumentation 
for the BBN test is shown in Figure 14. The interior microphone 
arrangements are shown in Figure 4. The BBN locations are 
not correct for measuring a true space average. The LaRC 
locations are correct in that sampling of equal volumes Is 
achieved. Unfortunately the BBN microphone axial locations were 
not discovered to be in error until after the tests were con- 
cluded. Comparison of BBN data to the LaRC data indicates that 
the erroneous axial locations did not introduce too much 
bias in the results (as will be seen shortly). 

For both LaRC and BBN tests, the cylinder was supported at its 
two ends via two soft suspension spring/cords. At BBN the 
experiment was conducted in a room with a volume of 420 m 3 
(14,820 ft 3 ) with dimensions 6 x 10 x 7 meters. The cylinder 
was approximately in the floor center about 1.5 meters above the 
floor. Scattering panels were positioned at the edges of the 
room and four JBL 4311 loud speakers, directed away from the 
cylinder, were used to develop a reasonably diffuse excitation 
field. The space-average excitation levels in the room were 
measured with a B&K 4134 condenser microphone, connected through 
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FIGURE 14. INSTRUMENTATION DIAGRAM FOR CYLINDER NOISE 
REDUCTION MEASUREMENTS WITH REVERBERANT 
FIELD EXCITATION 
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a B&K 2203 Sound Level Meter to a Spectral Dynamics SD-312 One- 
Third Octave Band Analyzer. The space-average levels were 
measured by traversing the microphone, mounted on a 2 meter boom, 
slowly over a path through the room volume. The interior sound 
levels were measured with similar instrumentation at radial and 
axial positions as indicated in Figure 4. Circumferential 
averages over four positions 90° apart were performed at each 
axial station using the "pause" facility available on the 
SD 312. Measurements were concentrated over one half of the 
cylinder since by symmetry one expects the interior level to be 
reflected about the center. 

Noise reductions, NFL , defined as the difference between the mean 
exterior sound pressure level and the interior level, at measure- 
ment position, i, were calculated and averaged in accordance with 
the relation 


NR = -10 log 10 


N 

1 V (-NR./10) 
N i 10 1 

i = l 


y 


where N is the number of sampled interior points. The result is 
the cylinder noise reduction, NR. The associated standard 
deviation for the noise reduction measurements was also calculated 
with 


s 


1 

N-l 


£ 

-i =1 


(10- N V 10 -10- NR/10) 


h 
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The BBN results are shown In Figure 15, together with ±ls 
limits. The lower bound lies closer to the mean curve than does 
the upper reflecting the energy averaging procedure used to 
calculate the mean, i.e., the higher interior levels (lower noise 
reductions) have greater contribution to the mean than do the 
lower levels (higher noise reductions). 

Figure 16 shows the LaRC measured noise reduction. Comparison of 
Figures 15 and 16 indicates rather consistent results. It is 
felt that the LaRC 4000 and 5000 Hz results are in error. At 
these frequencies, the LaRC noise reduction increases significant- 
ly. The two frequencies were the last two bands analyzed at 
LaRC. The BBN test specifically looked at the high end to 
determine whether this effect was real and although at 4000 and 
5000 Hz, the noise reduction increased, it remained below 10 dB 
as can be seen in Figure 15. Also at the extreme low frequency, 
(i.e., 40 Hz) a noise reduction of 50 dB was achieved through care- 
ful sealing of the cylinder in the BBN test. The noise reduction 
was measured to be 6 dB lower at LaRC. Nevertheless, the noise 
reduction curves are in good agreement. Both exhibit a trough 
at 125 and 160 Hz, a peak at 200 Hz, and then a general decrease 
to below 10 dB by 630 Hz. 

A more complete statistical analysis of the measured noise 
reductions is considered in Section 5, which deals with verifica- 
tion of the analytical model. 

4.3 Plane-Wave Tonal Excitation 


Sound levels were measured within the bare cylinder when excited 
by harmonic plane acoustic waves incident at angles of <p = 15° 
and 45° relative to the cylinder axis. The experimental set-up 


- 98 - 




FIGURE 15. BBN MEASURED NOISE REDUCTION (BARE CYLINDER) 
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was similar to that for the noise reduction tests except that the 
room was lined with absorptive foam, the cylinder was shrouded 
in foam wedges to prevent the development of a reverberant field 
in the room, and only one loudspeaker source, located about 
4 meters from the cylinder center, was used as the excitation. 
Figure 17 shows the test layout. A progressive wave field was 
produced along the cylinder axis. The equipment used for the 
tone transmission measurements is shown in Figure 18. An exterior 
microphone was located at the surface of the cylinder to monitor 
the excitation level, which was maintained constant to within 
±1 dB and equal to 90 dB (re 20 yPa) using the compressor circuit 
available on the B&K 1024 Sine Random Generator. The interior 
sound level was measured at 12 positions as shown in Figure 19, 
the microphone output being plotted on a B&K 2305 Graphic Level 
Recorder. Representative measured data for 15° and 45° angles 
of incidence are shown in Figures 20(a) and (b). All other 
measured data are presented in Appendix A . Table 4 summarizes 
the measurements which were made. 

It should be emphasized that the curves in Figures 20 and in 
Appendix A , are ratios of incident to interior pressure at 
specific locations. It was simply impractical to measure the 
ratio of incident to space-average interior, since many micro- 
phones would have had to have been used simultaneously. However, 
the interior levels from point to point have similar character- 
istics. By selecting two extreme incidence angles <J> = 15° and 45°, 
it is obvious that the transmission characteristics of the 
cylinder can change dramatically. For instance, the 15° and 45° 
cases look pretty similar out to about 400 to 500 Hz. However 
it is clear that the maximum interior levels (minimum ratio of 

<p? >/<p?>) occurs first at about 3000 Hz for 15° but much 

me i 0 

sooner, between about 1000 and 1600 Hz for 45 . At 2000 Hz the 


- 101 - 



Cylinder 



FIGURE 17. TEST LAYOUT FOR PLANE WAVE EXCITATION 
OF CYLINDER 
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FIGURE 18. EQUIPMENT SCHEMATIC FOR TONAL PLANE 
WAVE EXCITATION 
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FIGURE 19 . MEASUREMENT LOCATIONS FOR TONAL 
PLANE-WAVE EXCITATION 



105 


FIGURE 20(a) RATIO OF INCIDENT PRESSURE TO INTERIOR PRESSURE 
TONAL PLANE WAVE EXCITATION, 0 = 15 . 
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ratio of incident to interior pressure for the cases shown in 
Figure 20 is about 20 dB for 15° angle of incidence, and 
on its way down, while at 45° incidence angle, the ratio of 
incident to interior has reached its minimum and is on its way 
up. Thus the two selected incidence angles are sufficiently 
different to challenge the analytical model to duplicate the 
observed behavior. 


Table 4. Details of Interior Microphone Locations 
for Tone Transmission Measurements 


Rec. No. 

Axial 

Position 

(mm) 

Angle 

(degrees) 

Radial 

Position 

(mm) 

02 

203 

180 

208 

03 

203 

180 

127 

04 

203 

0 

127 

05 

203 

0 

208 

06 

406 

0 

208 

07 

406 

0 

127 

08 

406 

180 

127 

09 

406 

180 

208 

10 

610 

180 

208 

11 

610 

180 

127 

12 

610 

0 

127 

13 

610 

0 

208 

16 

610 

0 

208 

17 

610 

0 

127 

18 

406 

0 

127 

19 

406 

0 

208 

20 

203 

0 

208 

21 

203 

0 

127 


Angle of Incidence, 
<P = 15° 


Angle of Incidence, 
<f> = 45° 
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4.4 Point Mechanical Excitation 


Transfer functions between cylinder acceleration response and 
force input and between cylinder interior sound level and force 
input were measured on the test cylinder. Measurement procedures, 
instrumentation, and results are described below. 


The cylinder "inertance" will be defined as the difference 
between acceleration output and force input, viz.. 


so that 


H(u>) 


a(io) _ |a|e” la _ |a| -i(a-B) 

^ " | F |e- iB " l P l 

l H (“) I = jfrj- > 


is the inertance magnitude, with units of (mass”' 1 '), and 


<j>(w) = (a-B) 


is the inertance phase 
input locations may be 
"transfer inertance". 
function is defined by 

V“> 


angle. Acceleration measurement and force 
different; H(w) is then considered a 
The sound level-to-force input transfer 

= P( M ) _ M e~ lY jpj_ -i(Y-S) 

fTwT |F| -16 | F j 


so that 

|H s (w)| "Jfj-. > 

is the magnitude of sound-to-force transfer function, with units 
of (area -1 ), and N . 


<{> c (w) = (y-<S) , 

o 


is the corresponding phase angle. 
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The experimental configuration for these measurements is shown 
in Figure 21, and an instrumentation schematic is presented in 
Figure 22. Broad-band noise from the IE-20 Noise Generator was 
amplified and fed to the input of a Goodman's V50 Shaker. The 
shaker drove the test cylinder through a Wilcoxon Z602 Impedance 
Head which was glued with dental cement through a point stud 
attachment to the cylinder. The axis of the shaker was perpen- 
dicular to the cylinder skin. Figure 21 shows the locations of 
the single shaker force input (S), and of the five response 
measurement accelerometers (A1 to A5) and of the interior 
microphones (PI to P4). 

The impedance head contains its own (small) internal mass which 
loads the force gauge to produce a voltage even when the imped- 
ance head is unloaded (i.e., not attached to a test structure) 
and which is comparable to the force required to drive the test 
cylinder at, for example, cylinder resonance frequencies. 
Cancellation of this force due to the internal mass is required 
to give accurate transfer function measurements. A mass can- 
cellation circuit was inserted into the impedance head circuit 
to subtract from the force gauge output , a voltage propor- 
tional to the acceleration of the force gauge in such a manner 
as to minimize the force gauge signal when the impedance head 
was driven unattached to the cylinder. In this circuit, the 
accelerometer signal from the impedance head was inverted, 
proportioned, and subtracted from the force gauge output to 
minimize output signal. Figure 23 shows the signal spectra from 
the output of the mass cancellation circuit with and without the 
cancellation circuit in operation (driving only the impedance 
head); a 30 dB reduction in force output has been produced at 
most frequencies. 
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EXPERIMENTAL CONFIGURATION FOR POINT MECHANICAL 
EXCITATION OF TEST CYLINDER, SHOWING SHAKER AND 
MEASUREMENT LOCATIONS 











a 




Calibration of the force gauge was achieved by adding a small 
dummy mass (37- 5g) to the impedance head, driving at a known 
acceleration (lg) and noting the force signal level. The force 
and acceleration signals were fed to the Spectral Dynamics 
SD-360 Dual Channel Analyzer, operating in transfer function 
mode and the magnitude, phase, and coherence of H(w) were deter- 
mined and set for calibration. The sound level-to-force 
transfer function calibration was achieved by feeding the same 
voltage simultaneously to microphone and force sound level 
meters and adjusting the output levels. 

Representative transfer functions are presented in Figures 24 
and 25 respectively. All other measured data are presented in 
Appendix B. Table 5 summarizes the measurement data by record 
number, measurement, frequency range, and sensor location; 
reference should be made to Figure 21. 

Figure 24 shows the inertance magnitude | H ( f ) | 2 and the phase 
angle <p between the acceleration and force. Also shown is the 
coherence function y 2 which provides an estimate of the errors 
associated with the transfer function measurement . When y 2 ^1.0, 
the errors in | H | 2 and <t> are small but the errors increase as 
y 2 tends to zero. In general in the present measurements, 
y 2 > 0.2 and the estimates for | H | 2 are accurate to within 
about ±1 dB. 

Minima in the coherence are associated with maxima and minima in 
the acceleration response and to a lesser degree to spectral 
variations in the force input. At low frequencies the inertance 
magnitude | H | 2 increases with the cube of the frequency as 
expected for a stiffness-controlled response until the frequency 
of the lowest structural resonance. At this frequency the | H | 2 
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FIGURE 24. INERTANCE FUNCTION, H(f), FOR POINT MECHANICAL EXCITATION 
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Table 5. Summary of Measurement Information for Point 
Mechanical Excitation. 


Record 

No. 

Transfer 

Function 

Sensor 

Location 

Frequency 

Range 

102 



A1 

0 to 1 kHz 

103 



A1 

0 to 500 Hz 

104 



A2 

0 to 500 Hz 

105 



A3 

0 to 500 Hz 

105A 

Inert 

:ance 

A4 

0 to 500 Hz 

106 

I H ( f ) | 2 

A5 

0 to 500 Hz 

106R 

and y 2 

A5 

0 to 500 Hz 

107 



A5 

0 to 2 kHz 

108 



A4 

0 to 2 kHz 

109 



A2 

0 to 2 kHz 

110 



A3 

0 to 2 kHz 

111 

1 


A1 

0 to 2 kHz 

112 

j 


P2 

0 to 500 Hz 

113 



PI 

0 to 500 Hz 

114 

Sound level- 

PI 

0 to 2 kHz 

115 

to-f orce 

P2 

0 to 2 kHz 

116 

|H s (f)| 2 

P4 

0 50 2 kHz 

117 

<t> and y 2 

P3 

0 to 2 kHz 

118 



P3 

0 to 500 Hz 

119 

i 


P4 

0 to 500 Hz 
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spectrum peaks and the phase angle changes from l80° to 0°, 
both phenomena being quite consistent with the behavior of a 
single degree-of-freedom system at resonance. Above the funda- 
mental structure mode |h| 2 decreases in accordance with mass- 
controlled behavior, to a minimum value which coincidentally is 
the frequency of the lowest acoustic mode of the cavity. 

Successive peaks in | H | 2 and phase changes of ±180° are associ- 
ated with the presence of other structure modes. Comparison of 
the phase angle at different positions over the cylinder surface 
at particular frequencies reflects the modal structure of the 
acceleration response although such discrimination is only 
possible for lower frequency modes. Further it appears that all 
structure modes at most frequencies are driven by the point 
excitation and that circumferentially the antinodes of the modes 
lie along the axial surface line containing the excitation point. 

Figure 25 shows the magnitude and phase of the sound level-to-force 

transfer function, |H C f ) | 2 and <J> , respectively, as well as the 

s s 

coherence function y 2 for the measurement. As indicated, 

spectral peaks in |H | 2 correspond to the presence of acoustic 

s 

and structural resonant modes. Phase changes of 180° occur at 

each of these peaks and also at minima in | H | 2 , corresponding 

s 

to changes from mass to stiffness controlled cylinder response. 

The highest peaks in |H | 2 in a frequency band encompassing both 
structural and acoustic modes occur at the resonance frequencies 
of the acoustic modes, as, for example, at 142 Hz and 284 Hz, 
the (1,0,1) and (2,0,1) acoustic modes respectively. It is clear 
that the resonant structure modes do not couple as effectively 
to the interior acoustic field when poor circumferential match- 
ing in wavelength exists. 
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4 . 5 Progressive Wave Excitation 


Measurement of the sound transmission for the case of a random 
progressive wave traveling along the length of the test cylinder 
was the final test. The experimental configuration is shown in 
Figure 26. Four JBL 4311 Loudspeakers were placed in a plane 
perpendicular to the cylinder axis and approximately 3.5 meters 
from the center of the cylinder. These speakers were placed 
about 1 meter from the cylinder axis and at 90° positions. The 
speakers' axes were directed to point at the cylinder mid-length 
position; the effective angle-of-incidence was about 15°. The 
cylinder was supported in a position similar to that used in the 
previous plane wave tests and surrounded by sheets of absorptive 
foam on all sides to simulate an anechoic environment. Noise 
levels were measured with B&K 4133 12.5 mm condenser microphones 
close to the cylinder surface at four equi-spaced circumferen- 
tial locations (El to E4 ) at each of the three axial stations 
(z/L = 0.25, 0.50 and 0.75), as shown in Figure 26. Noise levels 
were measured inside the cylinder at four axial stations 
(z/L = 0.125, 0.25, 0.375, 0.50) with two interior microphones 
(B&K 4131 25.4 mm Condenser Microphones), at 127 mm (PI) and 
208mm (PO) radial distances. Circumferential averages over the 
0° , 90°, l80°, and 270° positions were carried out by use of the 
'pause' facility on the SD312 1/3 Octave Band Analyzer. Figure 27 
provides an instrumentation schematic. 

Initially the gains of the power amplifiers for each speaker 
were adjusted so that the sound level spectra at the circumfer- 
ential exterior locations (El to E4 ) at z/L = 0.50 were as 
similar as practical. Table 6 and Figure 28 contain sound level 
spectra measured in one-third octave bands which illustrate the 
variation in excitation levels, both circumferentially at each 
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Table 6. Exterior Sound Levels for Progressive Wave Excitation 


Axial Station (z/L = 0.25): 


Frequency 

Microphone El 
SPL 

E2 

SPL 

E3 

SPL 

E4 

SPL 

Mean 

SPL 

63 

89.50 

78.75 

87.20 

90.00 

89.00 

80 

90.75 

91.50 

87.50 

92.75 

91.00 

100 

92.75 

92.75 

89.75 

92.75 

92.20 

125 

89.00 

90.75 

89.50 

90.50 

90.00 

160 

84.50 

88.00 

83.OO 

85-75 

85.70 

200 

87.25 

74.50 

86.00 

89.00 

87.00 

250 

90.25 

78.50 

91.00 

80.50 

89.90 

315 

93-75 

79.50 

90.00 

79.50 

91.10 

400 

92.75 

80.50 

90.00 

81.50 

91.30 

500 

90.00 

83.25 

87.75 

83.00 

91.50 

630 

90.25 

80.50 

88.50 

81.50 

90.30 

800 

89.50 

80.50 

84.75 

81.50 

O 

O- 

-=r 

00 

1000 

88.75 

77.25 

89.75 

83.00 

90.20 

1250 

88.00 

76.50 

90.75 

81.50 

81.60 

1600 

O 

Lf\ 

OO 

CO 

79.75 

92.00 

81.50 

90.70 

2000 

89.25 

76.75 

91.00 

81.00 

89.80 

2500 

83.75 

74.50 

89.25 

86.25 

86 . 50 

3150 

83.OO 

81.50 

85.25 

86.50 

84.50 

4000 

82.25 

86.00 

86.75 

86.25 

85.60 

5000 

82.25 

85.50 

90.50 

90.00 

85.20 

6300 

84.75 

84.50 

88.25 

90.25 

87.60 

8000 

86.25 

82.50 

81.00 

89.25 

85.90 

10000 

83.OO 

81.00 

86.50 

87.00 

85.00 


- 121 - 





Table 6 . Cont ' d . 


Exterior Sound Levels for Progressive Wave Excitation 


Axial Station (z/L = 0.50): 
Microphone El 

E2 

E3 

E4 

Mean 

Frequency 

SPL 

SPL 

SPL 

SPL 

SPL 

63 

88.75 

88.00 

88.00 

87.75 

88.10 

80 

89.00 

89.OO 

85.50 

91.00 

89.00 

100 

92.25 

93-00 

88.00 

91.50 

91.60 

125 

89.00 

90.00 

88.75 

90.00 

89.50 

160 

86.00 

88.75 

83.50 

87.25 

86.80 

200 

86.00 

86.25 

82.25 

87.00 

85.70 

250 

88.25 

83.50 

88.25 

89.25 

87.80 

315 

90.50 

86.75 

87.75 

89.50 

88.90 

400 

92.25 

83.25 

87.50 

89.00 

91.10 

500 

91.50 

88.50 

86.75 

92.75 

90.50 

630 

86.25 

87.25 

85.75 

92.75 

89.00 

800 

88.00 

87.50 

85.00 

88.75 

87.50 

1000 

89.50 

88.75 

88.50 

91.75 

89.80 

1250 

87.25 

85.25 

89.25 

90.25 

83.50 

1600 

86.50 

89.50 

90.50 

91.50 

89.90 

2000 

88.25 

87.OO 

90.50 

90.75 

89-40 

2500 

83.50 

84.50 

89.25 

86.75 

86.60 

3150 

81.75 

82.75 

86.50 

85.OO 

84.40 

4000 

80.00 

85.OO 

87.OO 

87.75 

85.80 

5000 

80.75 

83.75 

89.50 

88.50 

86.90 

6300 

82.75 

84.00 

84.25 

90.00 

86.30 

8000 

84.00 

73.75 

84.75 

88.25 

85.60 

10000 

81.00 

72.75 

86.25 

88.25 

85.50 
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Table 6. Cont'd. 


Exterior Sound Levels for Progressive Wave Excitation 


Axial Station (z/L = 0.75): 
Micropnone El 

E 2 

E3 

E4 

Mean 

Frequency 

SPL 

SPL 

SPL 

SPL 

SPL 

63 

89.75 

89.75 

90.00 

90.50 

90.00 

80 

87.75 

87.50 

85.50 

90.00 

88.00 

100 

92.25 

92.75 

88.50 

89.50 

91.10 

125 

88.00 

89.50 

87.25 

85.75 

87.80 

160 

86.25 

86.50 

83.50 

87.00 

86.00 

200 

83.00 

88.50 

82.00 

86.50 

85.80 

250 

85.00 

85.50 

84.50 

87.75 

85.90 

315 

90.25 

78.75 

86.75 

87.25 

87.30 

-400 

91.75 

88 . 25 

86.75 

89.75 

89.50 

500 

89.75 

89.50 

83.25 

88.00 

88.30 

630 

85.75 

88.25 

83.75 

92.00 

88.60 

800 

88.00 

85.75 

83.75 

91.50 

88.20 

1000 

OO 

OO 

-O 

U~! 

85.50 

87.50 

90.75 

88.50 

1250 

87.00 

85.25 

87.25 

90.50 

87.90 

1600 

87.OO 

87.75 

89.25 

90.50 

88.80 

2000 

87.25 

87.25 

88.50 

91.00 

88.80 

2500 

82.75 

86 . 00 

87.75 

88.00 

86.60 

3150 

82.75 

82.00 

85.75 

85.50 

84 . 30 

4000 

83.50 

83.50 

88.50 

84.50 

85.60 

5000 

81.00 

86.25 

88.50 

87.75 

86 . 70 

6300 

80.75 

84 . 50 

86.50 

88.50 

85.80 

8000 

78.00 

84.00 

85.OO 

87.50 

84.70 

10000 

75.25 

82.75 

80.50 

87 . 00 

83.30 
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One-Third Octave Band Sound Pressure Level 
dB re 20/*.Pa 






axial station, and also between axial stations. The standard 
deviation between axial stations is generally less than 1 dB 
over the whole frequency range, while the standard deviation 
circumferentially is only marginally greater. The decay in 
axial level over the cylinder length is generally less than 1 dB 
based on the mean circumferential sound level spectra. Con- 
sidered from the viewpoint of excitation levels, the plane of 
speakers provides a circumferentially homogeneous excitation, 
with minimal axial decay in level. 

Coherence and phase measurements of the exterior pressure field 
were made over the cylinder surface, both circumferentially at 
z/L = 0.5 and along four equi-spaced axial lines. 

Figure 29 shows the locations of these measurement planes, while 
details are summarized in Table 7. The axial microphone separa- 
tions were 150mm, 300 mm, and 450mm, with the microphone closest 
to the speaker plane located at z/L = 0.25. The circumferential 
separations varied from 22^° (100mm) to 90° (400mm), and measure- 
ments concentrated on the top surface of the cylinder. Repre- 
sentative coherence and phase measurements are shown for axial and 
circumferential separations in Figures 30 and 31 respectively. 

For the axial coherence measurement (Figure 30), the phase varia- 
tion with frequency is essentially continuous with constant slope, 
characteristic of a non-dispersive traveling wave. The coherence 
is fairly high in a broad-band sense but varies strongly from 
frequency to frequency, reflecting contamination of the cross- 
spectrum measurement by multi-path transmission between trans- 
ducers, due, for example, to weak reflections from the absorptive 
surfaces . 
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4 Longitudinal Measurement Arrays 



FIGURE 29. MICROPHONE LOCATIONS FOR MEASUREMENT OF COHERENCE 
AND PHASE OF PROGRESSIVE WAVE FIELD 



Table 7. Summary of Coherence and Phase Measurements 


Run Figure Measurement Microphone Separation 

No. No. Direction Locations (mm) 


201 



(1,1) 

and 

(1,2) 



150 


202 



(1,1) 

and 

(1,3) 



300 


203 



(1,1) 

and 

(1,4) 



450 


204 



(7,1) 

and 

(7,2) 

Longitudinal 

150 


205 

30 

(7,1) 

and 

(7,3) 

( axial ) 

300 


206 



(7,1) 

and 

(7,4) 

Measurement 

450 


207 



(8,1) 

and 

(8,2) 



150 


208 



(8,1) 

and 

(8,4) 



450 


209 



(4,1) 

and 

(4,2) 



150 


210 



(4,1) 

and 

(4,4) 



450 


211 

d 


(1,3) 

and 

(4,3) 

i 

t 

400 

(90°) 

212 



(2,3) 

and 

(4,3) 

1 

Circumferential 

200 

(45°) 

213 

31 

(3,3) 

and 

(4,3) 

Measurement 

100 

(22.5°) 

214 



(3,3) 

and 

(4,3) 



100 

(22.5°) 

215 



(6,3) 

and 

(4,3) 



200 

(45°) 

216 



(5,3) 

and 

(4,3) 

1 


100 

(22.5°) 
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PROGRESSIVE WAVE TEST. 



For progressive wave excitation, the cross power spectral 
density of the external pressure field in the axial direction 
is assumed to be of the form 


<Cp(«.) = exp exp 

where E, is the separation distance (meters), U z is the trace 
velocity in the axial direction (meters/sec), and a z is the decay 
parameter. The trace velocity U z may be estimated from the 
average slope of the phase curve (Figure 30), using 


U 

z 


360 £ 

Slope 


3 


where the "Slope" is the frequency in Hz divided by the phase in 
degrees. Similarly the decay parameters may be estimated from 
the coherence y 2 (Figure 30). Since axially it is assumed that 
Y 2 = exp[-u(/(a U )], for large values of a (-50), a straight 
line approximation can be fitted to the coherence function 
giving 


cL 


Z 


2ttE 

U 

z 


2500 

log e Y f 




Appendix C contains all of the measured data. From these, a mean 

value of the axial decay parameter a =56.3 was determined and 

z 

a trace velocity U = 360m/sec was found. 

z 

In the circumferential direction, the phase plots do not show 
evidence of a progressive wave field. The circumferential field 
was assumed to have the characteristics of diffuse excitation 
which gives a coherence function of the form ( sinkn )/kn , where 
k is the acoustic wavenumber and ri is the separation distance 
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between observer points. Figure 31 shows the measured coherence 
in the circumferential direction together with the diffuse field 
approximation. These show some agreement. However, it is clear 
that the analytical representation is weak. 

Following the measurements of the coherence and phase, the space 
average interior levels were obtained. Measured exterior and 
interior sound levels are presented in Table 8 along with calcu- 
lated mean levels and standard deviations. The mean exterior 
excitation levels are also shown, together with the resulting 
mean differences between blocked and interior levels. One stan- 
dard deviation limits are also given. Where no upper limit to 
the difference is shown, the variation in the measured data is 
too great to allow a prediction of the upper limiting value. 

The spectra are plotted in Figure 32. 
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Table 8. Data for Progressive Wave Field Excitation 
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5.0 COMPARISONS AND VALIDATION OF ANALYTICAL MODELS 


Basically, validation is obtained through comparison of computed 
and measured results. An analytical model can be said to be 
validated if predictions made with the model are in agreement 
with experimental results, i.e., to some required degree of 
accuracy . 

One type of validation involves a comparison of a point pre- 
diction against a measurement having an uncertainty associated 
with it. The "true" value of the quantity being measured is 
known only within certain limits for any given level of confi- 
dence. The noise reduction measurement for the cylinder is an 
example of this type. Suppose that the average exterior level 
is known precisely. The measured space average interior level 
has associated with it an uncertainty due to the finite number 
of sampling points. The noise reduction computed from the 
measurements is an estimate of the true (real) noise reduction 
(which can be obtained only by sampling an infinite number of 
interior points). 

The variance of the interior measurements for each band can be 
calculated using the unbiased variance estimator 


N 



i~l 


where N is the sample size, x^ is the sound level at ith location, 

and x is the average of sound levels at N locations. Given the 

variance estimate for the measured space average levels in the 

cylinder, the critical value of the anticipated discrepancy 

between the predicted, x , and measured averages, say in each 

P ^ 

one-third octave band due to statistical variability is [29] 
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(63) 



/IT m ; a/2 


where s is the variance of the interior measurements, N is the 
sample size of the interior measurements, and t /n is the a/2 
percentage point of the Student "t" distribution with m degrees 
of freedom. The number of degrees of freedom is m = N-l. Eq.(63) 
simply defines the (1-a) percent confidence intervals. The 
critical values of the differences between the measured and 
predicted space-average interior levels, or, equivalently, the 
noise reductions, are then predicted with Eq.(63) by selecting, 
say, a = 0.01 for a 1 % level of significance. 

The interpretation of the result is as follows. If the true 
space-average levels inside the cylinder were exactly equal to 
the analytical predictions then there would be a probability of 
only 1 % that an indicated discrepancy between the measured and 
predicted levels based upon the experimental data would exceed 
the associated critical value. Hence, when an indicated dis- 
crepancy does exceed its associated critical value, i.e., lies 
outside the confidence intervals, it is very unlikely that the 
true space-average level and prediction are in fact equal, and 
one would accept the conclusion that a discrepancy actually 
exists. The indicated discrepancy is considered statistically 
significant at the 1 % level. On the other hand, if an indicated 
discrepancy does not exceed its associated critical value, the 
evidence is considered insufficient to assert that an actual 
discrepancy exists. 
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Unfortunately, this type of analysis really only gives an indi- 
cation of the accuracy of a given prediction, one which is develop- 
ed not just with an analytical model, but with an analytical model 
that utilizes input data that may or may not be correct. Thus 
in the overall scheme of things, what is being validated is not 
just the analytical model, but the ability to parameterize the 
model (with a reasonable effort) as well. For this reason, in 
considering the noise reduction, a best and final prediction is 
made for comparison against the measurements, i.e., through mani- 
pulation of the input data, a best estimate within the bounds of 
the available theory and input data is made and the quality of 
this prediction is compared, from a statistical standpoint, 
against the measurements. 

5 . 1 Basics of Computer Program 

The analytical models of Section 3 were computerized to automate 
the calculations of the cylinder noise reduction and the various 
transfer functions as required for comparison with experimental 
results. An outline of the overall computational procedure, 
with details of how the various components of the analytical 
model are assembled together within the computer program, is 
presented below. 

The initial step in all computations is to calculate, using the 
given dimensions, the resonance frequencies of all structure 

modes of the cylinder and to classify them into frequency bands. 

For the high-frequency regime, band-averaged response quantities 
including joint acceptances, modal densities, radiation effi- 
ciencies, etc., are then calculated over the frequencies of 
concern using appropriate excitation models as required. These 
results are calculated in a separate preparatory program and 
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stored on disc. The initial function of the main program is to 
calculate and order the resonance frequencies of the cylindrical 
volume in frequency bands, and to combine this acoustic modal 
data, including acoustic loss factors, with previously calculated 
and stored structural data (inluding structural loss factors). 

The detailed calculation of the cylinder noise reduction or 
transfer function is carried out with a series of subroutines 
which are structured in a similar manner to allow separation of 
resonance, stiffness, and mass-controlled structure and acoustic 
response. In particular, for the noise reduction, in each one- 
third octave band, the five modal pairs (structure and acoustic) 
which contribute most to the overall interior response are 
identified and rank-ordered. In addition, the contribution of 
resonant and non-resonant acoustic modes and of resonant and 
non-resonant structure modes to the interior sound level is 
determined in each frequency band. 

5 . 2 Noise Reduction of the Test Article 

Prediction of the test cylinder noise reduction is now considered, 
i.e., the development of the prediction which is used in the 
statistical comparison. This is followed by the statistical 
analysis itself. 

5.2.1 CaZg^Zated Ee^uVts 

The initial calculation was for the case of an ideal 1.2192 meter 
(48 inch) cylinder enclosing a 1.2192 meter (48 inch) cylindrical 
cavity. The input data consisted of the measured acoustic and 
structural loss factors as detailed in Table 3 • The noise 
reductions for the one-third octave bands with center frequencies 
from 40 Hz to 1250 Hz inclusive were computed with Eq.(10). 
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Beginning with the one-third octave band centered at 1600 Hz, the 
noise reductions were computed with the high frequency relation, 
Eq. (17). The selected division between low and high frequency 
regimes means that for predictions using the low frequency model, 
there were 30 or fewer acoustic modes resonant in each one-third 
octave band (for the cylinder of concern, the 1250 iiz band nas 
28 resonant acoustic modes ) . It should be kept in mind that in 
addition to comparing the predictions against measurements, one 
can also examine the transition from low to high frequency models. 

Figure 33 shows the initial prediction (Curve No. 1) plotted 
against the BBN measurement from Figure 15- As can be seen the 
prediction has a shape similar to the measurement. The predicted 
noise reduction exceeds the measured out to 315 Hz, but is less 
than the measured beyond that frequency until about 1600 Hz. Note 
that a negative noise reduction is predicted at 1000 Hz and at 
1250 Hz. Also there appears to be a step change in predictions as 
one moves from the low frequency to the high frequency models 
(from 1250 Hz to 1600 Hz). 

In order to determine if the step change was an artifact of the 
analysis, the cylinder length was taken to be 1.1938 meters 
(47 inches), which is more in keeping with the actual construc- 
tion of the test article. The cavity length remained 1.2192 
meters (48 inches). The resulting prediction is Curve No. 2. 

With this new length, the prediction for the 1250 Hz band changed 
from -10.2 dB to +3.0 dB. The 1000 Hz band, which for the 48 inch 
case was -5.3 dB, dropped to -11.8 dB. Predictions for other 
bands are, in comparison to measurements, in some instances 
better, in others, worse. The transition from the 1250 Hz band 
to the 1600 Hz band is now reasonably good. 
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INITIAL NOISE REDUCTION PREDICTIONS WITH 
ANALYTICAL MODEL 





As has been stated, the computer program can output the five 
highest contributing modal pairs. A sample of this output for 
the 1000 Hz band for Curve 2 of Figure 33 is given below. 

Noise Reduction for Frequency 1000 Hz 


Acoustic 

Mode 

Structure 

Mode 

P ( I NT ) **2/ 

Q N 

S 

Freq . 

M 

N 

Freq . 

P ( EXT ) **2 

1 2 . 

1 

671.3 

2 

2 

892.0 

.1780E-01 

3 2 

1 

780.4 

2 

2 

892.0 

. 2555E-01 

4 3 

1 

1063.9 

3 

3 

951.4 

.6043E-01 

0 3 

1 

902.9 

3 

3 

951.4 

.9543E-01 

2 3 

1 

945.8 

3 

3 

951.4 

.1501E+02 


P ( INT) **2/ NOISE 

TOTAL P (EXT) **2 REDUCTION 

.1525E+02 -11.83 


Contributions from Acoustic Modes 
Below Band In Band Above Band 

.4346E-01 .1519E+02 . 196 OE-OI 


It is obvious from this data that the large negative noise reduc- 
tion predicted is due to a single modal pair. Both of the 
problem modes are resonant in the 1000 Hz band and they are closely 
spaced. The (2,3,1) acoustic mode is resonant at 945.8 Hz; the 
(3>3) structure mode is resonant at 951.4 Hz. The analytical model 
thus predicts a 5.6 Hz difference in resonance frequencies. 

The modal bandwidths of the two modes can be calculated with the 
data of Table 3* For the acoustic mode (2,3,1) the half power 
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bandwidth is 


B = n f = .0097(9^5.8) = 9.17 Hz 
n n n 

The structure mode (3,3) half power bandwidth is 

B = n f = .0018(951.**) = 1.71 Hz . 

These are therefore two closely coupled modes. If these two 
modes were in fact truly closely coupled in the test article, 
the noise reduction might be less than that presently measured 
(perhaps negative). The interpretation of the measurements is 
that this close coupling does not occur in the test article. The 
coupling can be excluded in the computer model by forcing the 
difference between resonance frequencies to be greater than the 
modal bandwidths . 

In addition to the problem above, the loss factor n r for the 
structure modes used in the analytical model is the sum of the 
measured dissipative loss factor and the computed band average 
radiation loss factor. Examination of the principal structure 
mode contributers , such as the (3,3) mode in the 1000 Hz band, 
shows that from 630 Hz to 1250 Hz, all are acoustically fast 
modes and should have a radiation efficiency of unity. The 
average radiation loss factors which are used are five to six dB 
below the radiation loss factors of the principal contributing 
modes; thus the value of n is being underestimated. 

Curve 3 of Figure 33 is, again, the case of the 47 inch cylinder 
with a 48 inch cavity, except the spacing between the closely 
coupled modes (such as (2,3,1) and (3,3)) has been increased so 
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that there is no modal overlap. From 630 Hz to 1250 Hz, account 
has been taken of the higher radiation damping due to the princi- 
pal contributing structure modes being acoustically fast. This 
new prediction shows how important the input data are. If two 
closely coupled modes actually occur at low frequencies (in the 
real situation) the noise reduction will depend almost solely on 
the damping present. 

Proceeding with the analysis of the noise reduction curve, it is 
seen in Curve 3 of Figure 33, that between 400 Hz and 630 Hz, more 
transmission is predicted than measured. Shown below is the 
computer output for the bands from 315 Hz to 800 Hz. 


Acoustic 
Q N 

Noise Reduction 

for Frequency 315 

Hz 

: Mode 
S Freq. 

Structure Mode 
M N Freq. 

P ( I NT ) **2/ 
P ( EXT ) **2 

3 

4 

1 

1218.3 

2 

2 

331.0 

.7875E-04 

1 

4 

1 

1151.5 

2 

4 

331.0 

. 2353E-03 

0 

1 

1 

395.7 

1 

1 

692.9 

.7253E-03 

2 

2 

1 

714.2 

1 

2 

289.0 

.4580E-02 

0 

2 

1 

656.4 

1 

2 

289 .O 

.2791E-01 






P(INT)**2/ 

NOISE 



Total 


P(EXT) **2 REDUCTION 






• 3367E-01 

14.73 



Contri buti ons 

from Acoustic Modes 

Below Band 


In 

Band 


Above Band 

. 7077E- 

-07 


.1874E-04 


. 3365E-01 
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Noise Reduction for Frequency 400 Hz 


Acoustic Mode Structure Mode P(INT)**2/ 

Q NS Freq. M N Freq. P(EXT)**2 


1 

1 

1 

395.7 

3 1 

2349.7 

1178E-03 

2 

1 

1 

485.6 

1 1 

692.9 

3197E-03 

0 

2 

1 

656.4 

1 2 

O 

CT\ 

CO 

CM 

5150E-03 

1 

1 

1 

420.0 

2 1 

1661.6 

7326E-02 

0 

1 

1 

395.7 

1 1 

692.9 

5887E+00 





P ( INT ) * *2/ NOISE 




Total 

P ( EXT ) * *2 REDUCTION 





.5972E+00 

2.24 



Contributions 

From Acoustic Modes 


Below Band 

In Band 

Above Band 

. 1626E-06 

.5963E+00 

.9437E-03 


Noi se 

Reduction 

for 

Frequency 500 

Hz 

Acoustic Mode 
Q NS Freq. 

Structure Mode 
M N Freq 

P( I NT ) **2/ 
P ( EXT ) **2 

2 11 

485-6 

3 

1 

2349.7 

.1034E-02 

3 3 1 

996.7 

2 

3 

510.0 

. 16 09E-02 

Oil 

395.7 

1 

1 

692.9 

. 1779E-02 

1 — 1 
OO 

1 — 1 

913.8 

2 

3 

510.0 

.6124E-02 

2 11 

485-6 

1 

1 

692.9 

. 2184E+00 




P ( INT ) **2/ 

NOISE 

Total 


P(EXT) **2 REDUCTION 




. 2302E+00 

6.38 

Contr i buti ons 

from 

Acoustic Modes 


Below Band 


In Band 


Above Band 

.1811E-02 


. 2195E+00 

. 8953E-02 
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Noise 

Reduction for 

Frequency 630 

Hz 

Acoustic Mode 
Q NS Freq. 

Structure 
M N 

Mode 
Freq . 

P( I NT ) **2/ 
P ( EXT) **2 

2 11 

485.6 

1 

1 

692.9 

. 6262E-02 

Oil 

395.7 

1 

1 

692.9 

.1361E-01 

4 1 1 

688,0 

1 

1 

692.9 

.4765E-01 

12 1 

671.3 

2 

2 

892.0 

.9216E-01 

0 2 1 

656.4 

1 

2 

289.0 

.1075E+00 




P(INT)**2/ 

NOISE 

Total 


P(EXT)#*2 REDUCTION 




.2719E+00 

5.66 

Contr i buti ons 

from Acoustic Modes 


Below Band 


In Band 


Above Band 

.1988E-01 


. 2488E+00 


. 3169E-02 


Noi se 

Reducti on 

for 

Frequency 800 

Hz 

Acoustic 
Q N 

Mode 

S Freq. 

Structure Mode 
M N Freq. 

P ( I NT) **2/ 
P ( EXT) **2 

2 2 

1 

714.2 

3 

2 1491.4 

•3191E-02 

4 1 

1 

688.0 

1 

1 692.9 

. 3850E-02 

1 2 

1 

671.3 

2 

2 892.0 

. 7937E-02 

2 2 

1 

714.2 

1 

2 289.0 

.1908E-01 

3 2 

1 

780.4 

2 

2 892.0 

. 1642E+00 





P(INT) **2/ 

NOISE 


Total 


P(EXT)**2 REDUCTION 





.2124E+00 

6.73 


Contributions from Acoustic Modes 
Below Band In Band Above Band 

. 1416E-01 . 19HE+00 .7161E-02 
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The one thing that stands out in the data is that the (1,1) 
structure mode is the predicted dominant transmitter at 400 Hz 
and at 500 Hz. It also contributes significantly at 630 Hz. 

The structure model used is that of a freely supported cylinder. 
The (1,1) mode is a beam mode, i.e., a translational ring mode. 
However, the test article is suspended on soft spring/cords. If 
there is any mode which would be difficult to excite, it would 
be the (1,1) mode. In fact, an argument can be made for excluding 
all N=1 structure modes since they represent translation rather 
than bending in the circumferential direction. 

If one excludes the (1,1) mode, the predicted noise reduction is 
as shown in Figure 34. The measurements are more or less bounded 
from 400 Hz to 630 Hz by the predictions made with and without the 
participation of the (1,1) structure mode. 

Consider now the frequency range 125 Hz to 250 Hz. The computer 
output is 



Noise 

Reduct i on 

for 

Frequency 125 

Hz 

Acoustic 

Mode 

Structure 

Mode 

P( I NT ) **2 / 

Q 

N 

S 

Fr eq . 

M 

N 

Freq . 

P ( E XT ) **2 

1 

0 

1 

140.7 

2 

0 

3437.2 

•3597E-05 

4 

4 

1 

1273.9 

1 

4 

126.6 

. 4927E-05 

0 

3 

1 

902.9 

1 

3 

153.4 

. 1205E-04 

2 

4 

1 

1177.0 

1 

4 

126.6 

.1758E-03 

0 

4 

1 

1142.9 

1 

4 

126.6 

.7998E-03 






P ( INT) **2/ 

NOISE 


Total 


P (EXT ) **2 
.1006E-02 


REDUCTION 

29.98 


Contributions for Acoustic Modes 
Below Band In Band Above Band 


0 . 


3755E-05 


.9977E-03 


- 145 - 




31.5 63 125 250 500 1000 2000 4000 8000 

One-Third Octave Band Center Frequencies in Hz 


FIGURE 34. COMPARISON OF PREDICTED AND MEASURED 
NOISE REDUCTIONS 
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Noise Reduction for Frequency 160 Hz 

Acoustic Mode Structure Mode P(INT)**2/ 

Q N S Freq. M N Freq. P(EXT)**2 


2 5 

1 

1407-3 

1 

5 159.7 

. 4868E-04 

1 0 

1 

140.7 

2 

0 3437.2 

.5613E-04 

0 5 

1 

1378.9 

1 

5 159.7 

. 2136E-03 

2 3 

1 

945-8 

1 

3 153.4 

. 5755E-03 

0 3 

1 

902.9 

1 

3 153.4 

. 2813E-02 





P ( INT) **2/ 

NOISE 


Total 


P(EXT) **2 

REDUCTION 





. 3739E-02 

24 . 27 

Contributions 

for Acoustic Modes 


Below 

Band 


In Band 

Above Band 

. 6831E' 

-04 

0 

. 


. 3680E-02 


Noise 

Reduction 

for 

Frequency 200 

Hz 

Acoustic Mode 
Q NS Freq. 

Structure 
M N 

Mode 
Freq . 

P( I NT ) **2/ 
P ( EXT ) **2 

2 2 1 

714.2 

1 

2 

289 .O 

. 2624E-02 

pH 

on 

C\i 

945.8 

1 

3 

153.4 

. 2714E-02 

1 — 1 

on 

0 

902 . 9 

1 

3 

153.4 

.1331E-04 

021 

656.4 

1 

2 

O 

G\ 

CO 

C\J 

.1536E-04 

Oil 

395-7 

1 

1 

692.9 

. 1614E-04 




P ( INT) **2/ 

NOISE 

Total 


P(EXT) **2 REDUCTION 




. 5561E-04 

42.55 

Contributions from Acoustic Modes 


Below Band 

In 

Band 



Above Band 

.2026E-06 

0 . 




.5541E-04 
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Noise Reduction for Frequency 250 Hz 

Acoustic Mode Structure Mode P(INT)**2/ 

Q N S Freq. M N Freq. P ( EXT ) **2 


2 

1 

1 

485.6 


1 

1 692.9 

. 5411E-05 

2 

2 

1 

714.2 


1 

2 289.0 

.5245E-04 

1 

5 

1 

1386.0 


2 

5 264.0 

.6209E-04 

0 

1 

1 

395.7 


1 

1 692.9 

.7710E-04 

0 

2 

1 

656.4 


1 

2 289.0 

.3157E-03 







P(INT) **2/ 

NOISE 



Total 



P(EXT)**2 

REDUCTION 







.5244E-03 

O 

OO 

C\J 

CO 


Contri buti ons 

from 

Acoustic Modes 


Bel ow 

Band 



In 

Band 

Above Band 

.1052E 

-05 



0. 



.5243E-03 


Prom these data, it can be seen that the dip in the noise reduc- 
tion at 125 Hz is due to the resonant response of the (1,4) 
structure mode at 126.6 Hz. Note that an acoustic mode of N 
order equal to 4, occurring at 1142.9 Hz, and responding in a 
non-resonant fashion in the 125 Hz band, is the principal contri- 
butor. This shows clearly that it is no simple matter to make 
general statements regarding transmission at low frequencies. 

The prediction at 125 Hz is pretty much dependent on the loss 
factor of the (1,4) structure mode. 

At 160 Hz, the (1,3) structure mode is resonant and the noise 
reduction is predicted to be below 25 dB. At 200 Hz, the 
measured and predicted noise reductions are near 40 dB (LaRC 
measurement). In this band, there are no resonant acoustic or 
structure modes. 
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At frequencies below 125 Hz, the cylinder volume is stiffness 
controlled and the noise reduction is predicted to asymptote 
to about 53 dB, based on the volume displacement afforded by the 
end caps. 


At frequencies in the range 1250 to 5000 Hz, the noise reduction 
is predicted and measured to be quite low. The BBN measurement 
shows values as low as 3 dB. This indicates that radiation 
damping is playing a significant role in the energy dissipation 
of the structure modes. In this case, the transmission is 
resonance controlled. The mass controlled transmission terms 
in Eq. (17) are negligible, i.e., and Eq . (17) goes to 

I n 


<p 2 > , 

e s , t 


<D 7 > 


i s , t 


as , 

t rA 



ext 

*Vad 


Since the absorption inside is small compared to the radiation 
damping, the above becomes 


<p 2> , 

e s , t 

<p?> . 

i s , t 


n 


n 


ext 

rad 


Further with p <ti , , 
r rad 5 


it is found that 


<p 2 > , 

e s , t 

<p? > , 

s , t 


2 


or 


NR - 3 dB 

Apparently the correct explanation for the observed transmission 
at high frequencies (beyond 1250 Hz), is that the cylinder damp- 
ing is radiation controlled. 
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In Figures 35 and 36 the various modal contributions are shown. 
The interior pressure is always due to either stiffness or 
resonance controlled acoustic response (Fig. 35). At high fre- 
quencies, resonant structure modes are transmitting to resonant 
acoustic modes (Figures 35 and 36). 

5 . 2.2 Statistical Anal]i8i8 o£ C\jlindeY_ Noise Reduction Data 

As discussed in Section 4, considerable data were collected for 
the experimental cylinder in the form of exterior and interior 
sound pressure level measurements from which noise reductions in 
one-third octave bands were calculated. The data were the 
following : 

1) . LaRC tests - 40 noise reduction measurements using 

fixed interior microphone locations at 

(a) four angular locations (0 = 0, 90, 180, and 270 degrees) 

(b) two radial locations (r = 0.1397 and 0.2077 meters) 

(c) five longitudinal locations (z/L = 0.56, 0.69, 0.75, 
0.85, and 0.95). 

2) . BBN tests - 32 noise reduction measurements using fixed 

Interior microphone locations at 

(a) four angular locations (9 = 0, 90, 180, and 270 degrees) 

(b) two radial locations (r = 0.1397 and 0.2077 meters) 

(c) four longitudinal locations (z/L = O. 56 , 0.69, 0.75, 
and 0.88). 

Some of the data from the BBN test were machine averaged over the 
four angular locations so individual noise reduction estimates 
at all angles are not available. 
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One-Third Octave Band Center Frequencies in Hz 


FIGURE 36. STRUCTURE RESPONSE CONTRIBUTIONS TO 
RESONANT ACOUSTIC RESPONSE ONLY 
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To obtain an estimate of the cylinder noise reduction, a collec- 
tion of noise reduction measurements representing approximately 
equal acoustic subvolumes inside the cylinder is required. The 
angular and radial locations of the measurements for both the 
LaRC and BBN experiments fall reasonably close to the centers 
of equal subvolumes along these axes. The longitudinal measure- 
ment locations of the BBN tests do not represent equal subvolumes 
but they are sufficiently spread along the longitudinal axis to 
provide an acceptable approximation over half the cylinder length 
Symmetry permits these measurements to be used to represent the 
entire cylinder length. Since some of the angular measurements 
were machine average, the data at various angular positions for 
all the radial and longitudinal locations for both the LaRC and 
BBN tests were averaged to arrive at noise reduction estimates 
representing similar subvolumes for both tests. This provided 
a total of ten noise reduction estimates (two radial locations 
at five longitudinal locations) for the LaRC tests, and eight 
noise reduction estimates (two radial locations at four longi- 
tudinal locations) for the BBN tests. The resulting noise 
reduction estimates are presented in Tables 9 and 10. 


The noise reduction of the cylinder was computed using the 18 
estimates from the LaRC and BBN tests, assumed to represent 
equivalent subvolumes, as follows: 


NR = 


-10 log 


iX 

•i -1 


10 


-NR ± /10 


3 


where NR^ represents each of the N = 18 noise reduction estimates 
To quantify the accuracy of the estimates, 99% confidence inter- 
vals for the noise reduction estimates were also calculated using 
the relationship 
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Table 10. Noise Reduction Estimates from BBN Experiments 
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where 


NR 


99 


= -10 log 


-NR/1° ± B 

/N m;0, 


005 


j 


s 


1 

N-l 


f (io- NR i / 10 -io- NR/10 ) a 

1=1 


V 

2 


9 


and QQp. is the 0.005 percentage point of the Student "t" 

variable with m = N - 1 = 17 degrees-of-freedom. It is assumed 
that the individual noise reduction estimates are statistically 
independent, which is believed to be an acceptable assumption. 

The resulting space averaged noise reduction calculations and 
their 99% confidence intervals are detailed in Table 11 and 
plotted in Figure 37- Note that in two frequency bands (the 1/3 
octave intervals centered at 250 and 315 Hz), it was necessary 
to define the lower confidence limit for a 95% interval because 
a 99% interval produced an infinite upper bound. 

Also shown in Figure 37 is the analytical prediction for the 
noise reduction of the cylinder. In the 400 to 630 Hz region 
where the (1,1) mode of the cylinder dominates, the analytical 
result represents an arithmetic average of the case where (a) 
the (1,1) mode fully participates, and (b) the (1,1) mode does 
not participate. Note that the shape of the measured and pre- 
dicted noise reduction curves are similar, but the measured 
noise reductions generally exceed the predictions by statisti- 
cally significant amounts at frequencies below 160 Hz. The 
results at 125 and 160 Hz show that the loss factors of the (1,4) 
and (1,3) structure modes have likely been overestimated. In 
the volume stiffness controlled region, 100 Hz and below, the 
noise reduction clearly does not achieve as high a level as 


- 156 - 



Table 11. Estimated Mean Value of Noise Reduction 
and 99% Confidence Intervals 


1/3 Octave 
band center 
frequency 

(Hz) 

Mean Value 
of Estimated 
Noise Reduction 

(dB) 

99% Confidence 

Lower Bound 
(dB) 

Interval 

Upper Bound 
(dB) 

40 

45.7 

42.7 

62.2 

50 

45.7 

42.9 

56.0 

63 

42.7 

41.1 

45.3 

80 

43.2 

41.5 

45.9 

100 

39.9 

37.6 

45.2 

125 

17.8 

16.1 

20.6 

160 

16.3 

14.4 

19.4 

200 

36.6 

34.3 

4 l .6 

250 

14.9 

11.8 

25 - 7 * 

315 

17.6 

14 . 3 

25 . 4 * 

400 

12.9 

11.6 

14 . 9 

500 

14.9 

13.2 

17.6 

630 

8.1 

6.2 

11.6 

800 

8.0 

6.2 

11.0 

1000 

8.2 

6.6 

10.7 

.1250 

6.3 

4.5 

9.6 

1600 

7.4 

6.2 

9.1 

2000 

7.2 

5.8 

9.1 

2500 

5.4 

3.9 

7.5 

3150 

4.3 

2.7 

6.7 

4000 

7.0 

4.9 

11.4 

5000 

10.4 

7.8 

17.6 


* 95 % Confidence Limits 
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31.5 63 125 250 500 1000 2000 4000 8000 

One-Third Octave Band Center Frequencies in Hz 


FIGURE 37. 


MEASURED AND PREDICTED NOISE REDUCTION 
FOR CY LINDER 
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predicted. There was suspicion that the lap joint was the cause 
of the measured low frequency noise reduction being lower than 
predicted. When the lap joint effect was added to the computer 
model, the difference in the computed results with and without 
lap joint was less than 0.1 dB. 

5 . 3 Tone Transmission Comparison 

The computer mechanization for the tone transmission study is 
Eqs.(36) and (37). These give the difference between the inci- 
dent pressure level and the space average interior pressure 
level at frequency w . The predicted interior space averages are 
to be compared against the two composite sets for the 15° and 45° 
angles of incidence which are given in Appendix A , representing 
differences between incident pressure and interior pressure at 
various points (as opposed to the space average). 

Since the results in Figure 20 are typical, they are repeated 
here as Figures 38(a) and 39(a). Figures 38(b) and 39(b) are the 
computed results. The computations were made every 20 Hz from 
100 Hz out to 800 Hz, then in 50 Hz increments to 1400 Hz, and 
in 100 Hz increments to 2200 Hz. The predictions are in fairly 
good agreement with the measured results, basically exhibiting 
the measured behavior as discussed in Section 4.3* 

Table 12 lists the most significant contributing acoustic and 
structure modes at the various frequencies for the 45° angle of 
incidence. The data may be compared against the most signifi- 
cant contributing modes found for the noise reduction (band-limited) 
calculation as given in Section 5-2.1. At high frequencies, say 
1250 Hz and above, it is clear that coupled acoustic and structure 
modes lying close to the excitation contribute most to the 
interior level. 
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Table 12. Most Significant Contributions to <p?>/<p? nc > 
Plane Wave Excitation of (<|> = 45°) 


Freq . 

Acoustic 

q.N.S 

Mode 
Freq . 

Structural Mode 
M,N Freq. 

Modal 
C 0 n t r i b . 

Total 

100 

0,1,1 

396 

1,1 

693 

. 17E-6 

. 44E-5 

140 

1,0,0 

140 

2,0 

3437 

. 0033 

. 0034 

200 

0,1,1 

396 

1,1 

693 

. 114E-4 

. 179E-4 

240 

1 — 1 
C\J 

0 

656 

1,2 

289 

. 467E-4 

. 931E-4 

280 

0,2,1 

656 

1,2 

289 

. 00410 

. 00504 

300 

0,2,1 

656 

1,2 

289 

.00315 

. 00385 

O 

^r 

on 

0,1,1 

396 

1,1 

693 

.00058 

.00098 

360 

0,1,1 

396 

1,1 

693 

.00149 

. 00181 

O 

CO 

on 

0,1,1 

396 

1,1 

693 

. 00829 

. 00861 

400 

0,1,1 

396 

1,1 

692 

. 0742 

. 0748 

420 

1,1,1 

420 

2,1 

1662 

.162 

.165 

440 

0,1,1 

396 

1,1 

693 

. 00083 

.00150 

460 

2,1,1 

486 

1,1 

693 

. 00050 

. 00119 

480 

2,1,1 

486 

1,1 

693 

. 0116 

.0134 

500 

1,3,1 

914 

2,3 

510 

.00308 

. 00532 

520 

1)3,1 

914 

2,3 

510 

.00305 

. 00427 

540 

1 — 1 

on 

rH 

914 

2,3 

510 

.00047 

. 00089 

560 

4,0,1 

563 

3,0 

3437 

.00027 

. 00095 

580 

3,1,1 

579 

2,1 

1662 

. 00586 

.00762 

600 

1,2,1 

671 

2,2 

892 

. 00040 

.00071 

620 

1,2,1 

671 

2,2 

892 

. 00075 

. 00112 

640 

1,2,1 

671 

2,2 

892 

. 00190 

. 00345 

660 

1,2,1 

671 

2,2 

892 

. 0140 

. 0196 

680 

1,2,1 

671 

2,2 

892 

.0157 

. 0207 

700 

2,2,1 

714 

3,2 

1491 

. 00194 

. 00703 
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Most 

Plane 

Table 12. Cont'd. 
Significant Contributions to 
Wave Excitation of (<P = 45°) 

< D ? > / < D ? > 

K i K i nc 


Freq . 

Acoustic Mode 
q , N , S Freq . 

Structural Mode 
M,N Freq. 


Modal 
Contr i b . 


Total 

720 

2,2,1 

714 

3,2 

1491 


. 00737 


. 00887 

740 

2,2,1 

714 

3,2 

1491 


. 00043 


. 00118 

760 

3,2,1 

780 

2,2 

892 


on 

-=r 

0 

0 

0 


. 00128 

780 

3,2,1 

780 

4,2 

1966 


.1256 


.2052 

800 

5,1,1 

807 

4,1 

2725 


.00048 


. 00152 

850 

2,3,1 

946 

3,3 

951 


. 00117 


. 00329 

900 

0,3,1 

903 

3,3 

951 


.1750 


.1968 

950 

2,3,1 

946 

3,3 

951 

5 

. 000 

5 

.035 

1000 

3,3,1 

997 

4,3 

1376 


.0297 


. 0431 

1050 

4,3,1 

1064 

5,3 

1740 


. 0142 


. 0293 

1100 

4,3,1 

1064 

5,3 

1740 


. 00473 


. 0236 

1150 

1,4,1 

1151 

6,4 

1602 


. 0134 


. 0505 

1200 

2,4,1 

1177 

5,4 

1308 


. 0224 


. 0514 

1250 

4,4,1 

1274 

5,4 

1308 


.0631 


. 0823 

1300 

4,4,1 

1274 

5,4 

1308 


.5863 


. 64l4 

1350 

5,4,1 

1342 

6,4 

1602 


. 0420 


. 0606 

1400 

2,5,1 

14 07 

7,5 

1518 


.0397 


. 0606 

1500 

4,5,1 

1489 

7,5 

1518 

1 

.650 

1 

.818 

1600 

8,4,1 

1604 

7,4 

1855 


.1200 


.2578 

1700 

7,5,1 

1694 

8,5 

1739 

5 

.355 

5 

.838 

1800 

9,5,1 

1872 

8,5 

1739 


. 0126 


. 0410 

1900 

7,6,1 

1889 

10,6 

1848 


.1400 


.2373 

2000 

i — 1 

VO 

CO 

1966 

11,6 

2011 


.0862 


.1984 

2100 

11,5,1 

2073 

10,5 

2105 

5 

. 080 

5 

.165 

2200 

12,5,1 

2180 

11,5 

2255 


• 313 


.358 
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5.4 Mechanical Excitation 


Figures 40 and 4l show the predicted cylinder response for mechan- 
ical excitation plotted against the measured response. The 
measurements are for two different locations. Records 103 and 
104 (see Table 5 and Figure 21 for interpretation). 

Figures 42 and 43 give the predicted space average interior 
pressure for unit force input plotted on typical measured data 
for three interior points (Records 112, 118, and 119). Again 
Table 5 and Figure 21 should be consulted for the measurement 
locations. The predictions are in fairly good agreement with 
the measured data. 

5.5 Excitation by a Random Progressive Wave Field 

Before considering the results for excitation by a random 
progressive wave field, a few more words are required concerning 
the experimental data. 

It was found in Section 4.5 that the standard deviation was very 
large for the one-third octave bands from 100 Hz to 200 Hz 
inclusive (refer to Figure 32). A comparison of predicted and 
measured results is not feasible over this range. Sampling of 
many more interior points would be required to define the true 
space average interior level. However, the progressive wave 
test was an experiment outside the scope of work which was under- 
taken mainly for the purpose of obtaining additional experience 
with the cylinder for future consideration. Therefore, a 
prediction was made in spite of the experimental shortcomings, 
and is shown in Figure 44. As stated in Section 4.5, the cir- 
cumferential wave field was hard to model, i.e., the analytical 
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POINT MECHANICAL EXCITATION 


0 4000 8000 
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FIGURE 44. PREDICTIONS VERSUS MEASUREMENTS, BAND-LIMITED 
RATIO OF EXTERIOR PRESSURE TO AVERAGE INTERIOR 
PRESSURE. PROGRESSIVE WAVE EXCITATION 





representation is weak and therefore the prediction was not 
expected to be as good as it might have been had this represen- 
tation had a more solid basis. 

5. 6 Noise Reduction for Increased Interior Absorption 

In addition to the noise reduction tests discussed in Section *1, 
another experiment performed at LaRC considered the effect of 
additional interior absorption. Figure 45 shows the measured 
data for various points in the cylinder when a 0.00635 meter 
(0.25 inch) thick layer of Soundcoat acoustical foam was laid 
against the inside wall of the cylinder. The data are insuffi- 
cient to compute a space-average interior level and the actual 
noise reduction. However the noise reduction will lie close to 
the lower curves. Also plotted in Figure 45 is a BBN prediction 
of the noise reduction. Data from a Soundcoat Company brochure 
were used to estimate the acoustic loss factors (i.e., absorp- 
tion coefficients) needed in the analysis. Since the damping of 
the cylinder structure modes was unknown, the prediction is based 
on a value n r = 0.01 for all bands, i.e., it was assumed that the 
addition of the foam would increase the structure loss factors 
to at least that level. A comparison is not intended other than, 
perhaps, the statement that in general the same trends are shown. 
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APPENDIX A 


Plane Wave Tonal Transmission Data 
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APPENDIX B 

Mechanical Excitation 


Data 
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APPENDIX C 

Progressive Wave Field Data 
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